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PREFACE 
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I. INTRODUCTION TO MARSYAS 




A. General Description of MARSYAS 

MARSYAS (Marshall System for Aerospace Simulation) was developed by 
NASA's Computation Laboratory at Marshall Space Flight Center to furnish 
engineers with a software system that allows quick and easy simulation 
of physical systems on a digital computer. 

MARSYAS is a simple, flexible language which can be coded by users 
who are unfamiliar with computer programming. It is designed for the 
engineer with little experience in simulation who desires to simulate 
large physical systems. The language can be used to solve a system of 
differential equations or to simulate control systems including analog 
computer block diagrams or both simultaneously. Thus, the user has the 
ability to mix differential equations with diagrams in his model. The 
block diagrams can contain, among other things, adders, integrators, 
transfer functions, multiple input/output nonlinear devices, algebraic 
equations and nonlinear ordinary differential equations. A block dia- 
gram is specified by the user-given names of its models and submodels, 
inputs and outputs, element names, parameters (if any), and their inter- 
connections. Submodels can be nested to any degree required. With 
MARSYAS, no preset pattern of connecting elements is required. Elements 
can be connected in pairs, groups or any manner desired by the user. A 
large library of Standard Elements and Excitation Functions is part of 
the MARSYAS system. DEVICE and FUNCTION statement operators allow the 
user to construct unusual element or excitation functions as needed. 
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I A. GENERAL DESCRIPTION OF MARSYAS 

MARSYAS is a flexible language in that, with few exceptions, there 
is no rigid statement operator structure within a given module. Most 
statements can be used without regard for the order in which they appear 
within the modules. Depending upon the computer system in which MARSYAS 
is installed, the user has the capability of storing models in a Func- 
tional Data Base. The Fortran Object Program generated from the MARSYAS 
source program can be extracted and run separately, if the user's com- 
puting facility can accommodate this feature. When using CHANGE opera- 
tors, the user has multiple simulation capability without the necessity 
of either rewriting his model or resubmitting his deck. 

An elaborate plotting system is part of the MARSYAS language allow- 
ing the user nearly unlimited flexibility in specifying his graphical 
output. Additionally, the Fast Fourier Transform of any output variable 
can easily be obtained. A tabular listing of a model in the Functional 
Data Base or of a model currently being run can be obtained using the 
LIST operator. 

Automatic features of MARSYAS include the detection and solution 
of linear and nonlinear algebraic loops. For problems which contain 
discontinuities, the MARSYAS system automatically changes integration 
schemes to integrate through the discontinuity, unless instructed other- 
wise by the user. 

MARSYAS is designed in modular form so that modifications to the 
system models can be made with a minimum of effort. In order to achieve 
comprehensive analysis capability and effective computation, modern con- 
trol theory is used as the mathematical foundation of MARSYAS. The 

8 . . 
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I A. GENERAL DESCRIPTION OF MARSYAS 

differential equations generated from block diagrams, or coded as equa- 
tions, are rearranged internally into vector-matrix state equations 
which are then solved . 

The language is designed so that the user transmits to the computer 
only the information essential to describe the mathematical model and 
specify the simulation run. 

MARSYAS is divided into four successive modules which describe inde- 
pendent functions of the simulation. These modules are as follows: 
Description Module 
Modification Module (optional) 

Simulation Module 
Post Processing Module 

The user has the ability to control some of the internal processing 
of the simulation by specifying his numerical integration method, inte- 
gration step size or even the truncation error. Normally, he need not 
concern himself with these details since MARSYAS handles these details 
automatically. 

MARSYAS names can be up to 36 characters in length so that the same 
names as found in engineering documentation can be used. The MARSYAS 
alphabet consists of the letters A through Z, the numbers 0 through 9, 
and the backward slash (\ ). There are no reserved words in MARSYAS. 
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I Bo THE MARS Y AS PROGRAMMING SYSTEM 


B. The MARSYAS Programming System 

The MARSYAS programming system consists of two basic components: 

1. A Source Program which consists of a set of MARSYAS statements 
which may contain a set of FORTRAN subprograms. 

2. A Processor Program which pre-compiles the MARSYAS language 
into a set of FORTRAN programs called the Object Program. 

The Source Program is fed to the computer on cards. The MARSYAS 
coding format is "free form" and information may be punched into any 
card column, using as many cards as necessary to complete a statement. 
Statement operators are always followed by a colon (:) and ended by a 
dollar sign ($). A colon used without a statement operator repeats the 
previous operator. Imbedded blanks in the coding are ignored. If FORTRAN 
subprograms are used, the coding must adhere to the rules of FORTRAN. 

The MARSYAS Processor Program converts the MARSYAS source language 
into FORTRAN code. This FORTRAN code, also called the Object Program, 
is then processed and executed by the computer operating system in the 
same manner as any other FORTRAN source coding. 

The MARSYAS processor program is written in FORTRAN V and, at the 

\ 

present time, can be implemented only on a Univac 1108 computer. Since 
the processor is modular in design, it can be modified so that the full 
system can be installed on other large computers. 
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I C. MARSYAS PROGRAM CONTROL STATEMENTS 

C. MARSYAS Program Control Statements 

There are two program control statements necessary to control the 
execution of the MARSYAS program. The BEGIN statement identifies the 
MARSYAS system model being executed and the END statement marks the end 
of the MARSYAS coding. The use of the BEGIN and END statements is illus- 
trated in Figure 1. The names given to the system model must be identical 
for the BEGIN and END statements. The general format for both is as 
follows : 

BEGIN: name $ 

MARSYAS Program 
END : name $ 

The LIST operator is a third program control statement used to cre- 
ate a tabular listing of a MARSYAS model currently being run or previously 
stored in a Functional Data Base (if any). 

The LIST statement will create a list of the following: 

a) The model name. 

b) Model input and output terminals with alternate names, if any. 

c) Element mnemonics, names and alternate names, if any. 

d) Parameter names and associated values. 

e) Connections between elements and system input-output terminals. 

f) A list of submodels with their inputs and outputs. 

The general format for the LIST operator is as follows: 

LIST: model name $ 

The LIST statement may be placed anywhere after the BEGIN statement 
and before the SIMULATE statement, even within a Description or Modifica- 
tion Module. 
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ID. THE STRUCTURE OF THE MARSYAS LANGUAGE 


D. The Structure of the MARSYAS Language 

Usually, engineers prefer to describe the system being simulated 
using block diagrams since this form of "graphical" representation is 
visually comprehensive. The blocks of the diagram can have multiple 
inputs and multiple outputs and blocks imbedded within blocks. The low- 
est level block is called an ELEMENT, the highest, a MODEL. A SUBMODEL 
is a model imbedded within a MODEL. 

A Description Module is used to describe the structure of a model 
given in block diagram or equation form. It is headed by the operator 
MODEL and terminated by an END$ statement. The ELEMENTS statement con- 
tains the name of the element, its mnemonic, and its parameters. 

ELEMENTS are devices which may be linear or nonlinear. A linear 
element can be as simple as a constant multiplier or as complex as a 
transfer function. Nonlinear elements are representations of either 
algebraic equations or switching functions or memory devices such as 
hysterisis. Frequently-used elements are listed in the Table of Standard 
Elements at the back of the Manual. If an element needed is not found 
in the Table of Standard Elements, a FORTRAN subroutine called a DEVICE 
can be constructed to form the needed element. Parameters are con- 
stants written in the format shown in the Table of Standard Elements 
and are either numerical values or names. The numerical value of a named 
parameter is given by the PARAMETERIZE Statement. The CONNECT State- 
ment connects strings of inputs and outputs of elements, submodels, 
system inputs, or system outputs, to form the system block diagram. For 

12 
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FIGURE 1 SIMULATION OF A MATHEMATICAL MODEL USING MARSYAS 
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ID. THE STRUCTURE OF THE MARSYAS LANGUAGE 

elements or submodels having a single input and output, only the name 
of the element or submodel appears in the CONNECT Statement. The INPUTS 
Statement designates names of the inputs of the model; the OUTPUTS State- 
ment designates names of the outputs of the model. 

The Modification Module allox^s inserting, deleting, and disconnecting 
of elements and submodels previously described in the Description Module 
through the use of the SUBSTITUTE, DELETE and DISCONNECT Statements. 

The Simulation Module completes the specifications of the system 
being analyzed. The INITIALIZE Statement specifies the initial condi- 
tions for the integrators and transfer functions. Excitation functions 
applied at the system inputs are specified with EXCITE Statements. Fre- 
quently used excitation functions are listed in the Table of Standard 
Excitation Functions at the back of the Manual. If an excitation func- 
tion needed is not found in this Table, a FORTRAN subprogram called a 
FUNCTION can be constructed to form the needed excitation. The INTEGRATE 
Statement specifies the integration method to be used if a method other 
than the standard method is desired. The STOP IF and TERMINATE 
IF Statements determine the condition (s) under which the simulation is 
halted. A CHANGE Statement is used for performing repetitive simulations 
without the need for dismantling the system model. 

The Post Processing Module follows the Simulation Module and speci- 
fies the format of data presentation to the MARSYAS processor. Both 
tabular data and graphs are available for presenting the output data by 
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I D. THE STRUCTURE OF THE MARSYAS LANGUAGE 

specifying the PRINT and PLOT statements, respectively. The printing 
interval is specified in the SAMPLE Statement. A- frequency analysis can 
be obtained at any output terminal through the use of the FOURIER State- 
ment. See Section II for a more detailed description of each of the 
Modules mentioned in this section. 

E. How to Use MARSYAS - Basic Ideas 

The complete simulation of a model under the MARSYAS system is com- 
posed of several phases or modules as shown in Figure 1. The following 
brief description of Figure 1 will enable the user to follow the con- 
struction of the simple example shown on the following pages. A detailed 
explanation of the modular structure of MARSYAS can be found in Section II. 

The BEGIN statement identifies the user-given MARSYAS program 
name (Example: BEGIN: TEST RUN1$). 

The optional FORTRAN DEVICE and FUNCTION statement box, which ap- 
pears at the top of Figure 1, is used for placing unusual excite and ele- 
ment descriptions not found in the standard MARSYAS library. In general, 
the user-coded FORTRAN DEVICE and FUNCTION subprograms may appear any- 
where within the MARSYAS program provided they appear before being 
referenced in the MARSYAS coding. The safest approach, however, is to 
place such subprograms ahead of the Description Module as shown to avoid 
difficulty. 

The Description Module is used to define the structure of the model. 
There may be more than one Description Module in a MARSYAS program. It 
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I E. HOW TO USE MARSYAS - BASIC IDEAS 


may be modified in a Modification Module. When using the block diagram 
problem formulation, six distinct statement operators must be specified, 
in any order with the exception of MODEL and END which must be first and 
last, respectively, as follows: MODEL name, INPUTS, OUTPUTS, ELEMENTS, 

CONNECT statements and END. The MODEL name is user-defined and is ref- 
erenced by the Simulation Module (Example: MODEL: MARSYAS ILLUSTRA- 

TION MODEL $) . The INPUTS and OUTPUTS statements define all the model 
inputs and outputs. If the model has one input and one output whose 
user-assigned names are U1 and Y1 , respectively, the statements would 
read as: 

INPUTS: Ul$ 

OUTPUTS: Y1 $ 

The ELEMENTS statement (s) specifies all of the elements contained 
within the model, such as adders, integrators, transfer functions, etc. 

A complete list of MARSYAS elements appears in the Table of Standard 
Elements, along with their mnemonics. For example, the mnemonic for an 
integrator is IN and the mnemonic for an adder is AD. If the model con- 
tains one adder and two integrators, the ELEMENTS statement would read: 

ELEMENTS: AD, ADDER1 $ : IN, INTEGRAT0R1, INTEGRAT0R2$ 

The CONNECT statement (s) connects the elements, inputs and outputs 
of the model to form a complete circuit. Using the inputs, outputs, and 
elements referred to above, the CONNECT statement could read: 

CONNECT: Ul, INTEGRATOR!, ADDER1 , INTEGRAT0R2 , Yl$ 
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IE. HOW TO USE MARSYAS - BASIC IDEAS 


The CONNECT statement is easy to visualize in the following sche- 
matic : 



INTEGRATOR! ADDER! INTEGRATOR2 


For the simple model just shown, the complete Description module 
is as follows: 

MODEL: MARSYAS ILLUSTRATION MODEL$ 

INPUTS : Ul$ 

OUTPUTS : Yl$ 

ELEMENTS: AD, ADDERl$ 

: IN, INTEGRATORl , INTEGRATOR2 $ 

CONNECT: Ul, INTEGRATORl, ADDER1 , INTEGRAT0R2 , Yl$ 

END$ 

The END statement is the last statement to appear in the Description 
Module. Every Module shown in Figure 1 must terminate with an END state- 
ment. 

The Simulation Module follows the Description Module. This module 
begins with the SIMULATE statement. It references the main model name 
being simulated as: 

SIMULATE: MARSYAS ILLUSTRATION MODEL$ 
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I E. HOW TO USE MARSYAS - BASIC IDEAS 




SIMULATE must be the first and END the last statements in the Simu- 
lation Module. The remaining statements described below may be located 
anywhere within the Module, without regard to order. The Simulation 
Module contains the following statements: EXCITE, INTEGRATE (optional), 

TERMINATE IF or STOP IF, INITIALIZE, CHANGE (optional) and END. 

The EXCITE statement is used to specify an excitation function on 
a particular input. Usually, but not always, a Standard MARSYAS excita- 
tion function is used. If an excitation function not listed in the Table 
of Standard Excitation Functions is needed, the user may construct a 
FORTRAN subprogram to do the job. This subprogram is called a FUNCTION 
and is placed before the Description Module as previously discussed. 

A typical EXCITE statement might read: 

EXCITE: Ul, FSIN (2.0, 3.0, 4.0)$ 

where FSIN is the MARSYAS mnemonic for a sinusoid. The above EXCITE 
statement says, "Excite input terminal Ul with 2 sin (3t + 4)." 

The INTEGRATE statement specifies the mode of integration to be used 
in solving the model. If the INTEGRATE statement is omitted, the problem 
will be solved using the Sarafyan variable-step method. If fourth-order 
Runga-Kutta is desired, the INTEGRATE statement could read: 

INTEGRATE: RK, TIMESTEP, 0.01$ 

where RK is an abbreviation for Runga-Kutta and TIMESTEP, 0.01 specifies 
the integration interval. If the word TIMESTEP is omitted, a step size 
of 0.01 is used automatically. 
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I E. HOW TO USE MARSYAS - BASIC IDEAS 


It is usually necessary to terminate the problem solution when a 
certain condition has been reached in one of the model parameters 
or when a certain point in time is reached. If, for example, it is 
desired to stop the solution after ten time units, one of the following 
statements could be used: 

TERMINATE IF: TIME .GT. 10.0$ 

or 

STOP IF: TIME .GT. 10.0$ 

The INITIALIZE statement is used to impose initial conditions, 
other than zero, on integrators and transfer functions. If, for example 
an initial condition of 3.2 were required on INTEGRAT0R1 in the above 
example, the statement would read 

INITIALIZE: MARSYAS ILLUSTRATION MODEL, INTEGRATOR! (3.2)$ 

When using the INITIALIZE statement, the model name as well as the 
element name must be specified. If there are no initial conditions asso 
ciated with the problem, the INITIALIZE statement is omitted. 

When it is necessary to make a succession of runs with different 
parameter elements, a CHANGE statement may be used. CHANGE statements 
are discussed in detail in Section II C. 
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I E. HOW TO USE MARSYAS - BASIC IDEAS 


The END statement is the last statement to appear in the Simulation 
Module. Summarizing the statements discussed above, the Simulation 
Module for the simple example is as follows: 

SIMULATE: MARSYAS ILLUSTRATION MODEL $ 

EXCITE: Ul, FSIN (2.0, 3.0, 4.0)$ 

INTEGRATE: RK, TIME STEP, 0.01$ 

TERMINATE IF: TIME .GT. 10.0$ 

INITIALIZE: MARSYAS ILLUSTRATION MODEL, INTEGRATOR1 (3.2)$ 

END$ 

The POST PROCESSING Module follows the Simulation Module and is used 
to specify which outputs to record, the time interval and the method of 
presentation for viewing, that is, tables and/or graphs. In addition, a 
frequency analysis can be specified at any output using the FOURIER state- 
ment. The following statements are used in the POST PROCESSING Module: 
PRINT, SAMPLE, PLOT, FOURIER, END. An example of the use of these state- 
ments is as follows: 

PRINT: Ul, Y1 $ ' 

SAMPLE: STEP, 10$ 

PLOT: LINEAR (0.0, 10.0, 1), Ul, Yl$ 

END$ 
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I E. HOW TO USE MARSYAS - BASIC IDEAS 

As a consequence of the above statements, the MARSYAS processor 
will tabulate Time, U1 , Y1 at every tenth integration step (STEP, 10) 
and plot on separate sheets of 1 ine ar graph paper (SC 4020 plotter) U1 
and Y1 versus time in the time interval 0-10. As usual, the END state- 
ment is the last to appear in the Module. 

The very last statement to appear in the MARSYAS coding is a final 
END name statement as: 

END: TEST RUN! $ 

This is the counterpart of the BEGIN statement which is the first 
statement to appear in the MARSYAS deck. Its purpose is to signal the 
MARSYAS processor that there is no more MARSYAS coding. 

Combining all of the above MARSYAS statements as shown in Figure 2 
will emphasize the basic structure of the language and will enable the 
user to follow the solution of the sample problem on the following pages. 
A detailed description of all of the statement operators discussed in 
this section will be found in Section II of the Manual. 
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PROCESSING , 

MODULE SIMULATION MODULE DESCRIPTION MODUL 


IE, NOW TO USE MARSYAS - BASIC IDEAS 

BEGIN: TEST RUN1$ 

K 

‘ MODEL: MARSYAS ILLUSTRATION MODELS 
INPUTS: U1$ 

OUTPUTS: Y1S 
ELEMENTS: AD, AODER1S 

ELEMENTS: IN, INTEGRATORS INTEGRATOR2S 
CONNECT: U1, INTEGRATOR!, ADDER!, INTEGRATOR2, Yl$ 

ENDS 

7 

* 

SIMULATE: MARSYAS ILLUSTRATION MODELS 
EXCITE: Ul, FSIN (2.0, 3.0, 4.0}$ 

INTEGRATE: RK, TIMESTEP, 0.0!$ 

TERMINATE IF: TIME .GT. 10.0$ 

INITIALIZE: MARSYAS ILLUSTRATION MODEL, INTEGRATOR! (3.2)$ 
ENDS 

F— 

PRINT: Ul, Yl$ 

SAMPLE: STEP, 10$ 

LINEAR (0.0, 10.0, 1), Ul, Y!$ 

END: TEST RUN1S 


FIGURE 2. ILLUSTRATION OF MARSYAS CODING 
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I F, A SAMPLE PROBLEM 


F. A Sample Problem 

Consider the simple mass-spring-dashpot system shown in the illus 
tration below,. 



The motion of the mass is described by the following differential 
equation and initial conditions: 

Mx + Bx + Kx = tJ (1) 

x(0) = 20.0 and x(0) » 0 (2) 

where M = 10.0, B = 2.5, K ~ 8.6, and U 3 sin t 

There are two methods of solving this problem using MARSYAS: 

(i) Equation (1) can be .transformed into a block diagram whose elements 
will consist of two integrators, an adder, and three constant multipliers 
or, (ii) the equation can be coded directly as shown above. The two 
methods of coding are illustrated in Figures 3 and 4, 

The EQUATION option appears to be the more straight-forward method 
of solution for this particular problem in that the Description Module 
requires less coding. However, larger systems would probably be solved 
using the block diagram approach since more insight into the structure 
of the problem is thereby gained. 
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IF. A SAMPLE PR03LEM 



BEGIN: MARSYAS EXAMPLES 
MODEL: MASS SPRING DAMPERS 
INPUTS: US 

OUT PUTS: AD TOUT, XDOT, X$ 

ELEMENTS: IN, INI, !N2$ : AD, AD1$ : CM, CM1{2.5), CM2(0.10), CM3(8.6)$ 


CONNECT: ADI, ADIOUTS : U, ADI, CM2, INI, IN2, CM3, -ADIS : INI, XDOTS 
: INI, CM1, -AD1$ : IN2, X$ 


! 
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LU 
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o 
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ENDS 

SIMULATE: MASS SPRING DAMPERS 
EXCITE: U, FSIN (1.0, 1.0, 0.0)$ 

INTEGRATE: RK, TIMESTEP, 0.01$ 

STOPIF: TIME .GT. 5.0$ 

INITIALIZE: MASS SPRING DAMPER, IN2(20.0)$ 
ENDS 

PRINT: X, XDOT, ADIOUT, US 
SAMPLE: STEP, 1$ 

ENDS 


END: MARSYAS EXAMPLES 


FIGURE 3. BLOCK DIAGRAM APPROACH TO SOLVING A SECOND 
ORDER PROBLEM USING MARSYAS 
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IF. A SAMPLE PROBLEM 


BEGIN: MARSYAS EXAMPLES 

5 

MODEL: MASS SPRING DAMPER, EQUATIONS 
INPUTS: U$ 

OUTPUTS: AD1QUT, XDOT, X\OUT$ 

EQUATION: 10.0* X“ + 2.5 * X‘ + 8.6 * X = US 
: ADIOUT = M * X”$ : XDOT - X‘$ 
s X\OUT = XS 

ENDS 

— 

SIMULATE: MASS SPRING DAMPERS 
EXCITE: U, FSIN (1.0, 1.0, 0.0)$ 

INTEGRATE: RK, T1MESTEP, 0.01$ 

STOPIF: TIME . GT. 5.0$ 

INITIALIZE: MASS SPRING DAMPER, X(20.0)$ 
ENDS 

| PRINT; X\OUT, XDOT, ADIOUT, U$ 

SAMPLE: STEP, 1$ 

| ENDS 

END: MARSYAS EXAMPLES 


FIGURE 4. EQUATION APPROACH TO SOLVING A SECOND ORDER 

PROBLEM USING MARSYAS 
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IF. A SAMPLE PROBLEM 


The program control statements BEGIN name and END name must appear 
regardless of which method of coding is selected.- The name MARSYAS 
EXAMPLE is selected to identify this MARSYAS simulation. 

The following explanation of the coding which appears in Figures 3 
and 4 should enable the user to become familiar with the MARSYAS language. 
A complete description of each of the Modules -and Statement operators 
appears in Section II of this Manual. 

Description Module 

The name MASS SPRING DAMPER identifies the MODEL whose structure is 
being described in the statements to follow. If SUBMODELS were referenced 
in this MODEL, then additional Description Modules would be required to 
describe the structure of each submodel. Notice that in Figure 4 the word 
EQUATION follows the model name indicating that equations are to be used 
in the MODEL. The excitation U is an input to the system and is so stated 
in the INPUTS statement. Since we wish to observe the behavior of the 
system being excited, the names X, XDOT and AD10UT have been specified as 
outputs. Generally, outputs can be at any point in the system so long as 
they are specified in the OUTPUTS statement. 

When the problem is simulated in block diagram form, all of the 
elements and connections, as well as the inputs and outputs, must be 
specified in a Description Module as shown in Figure 3. The constants 
B, M and K are constant multiplier elements with names CML, CM2 and CM3, 
respectively. The numerical value of each constant multiplier appears 
in parenthesis after its name in accordance with the format shown in 
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the Table of Standard Elements. The two integrators in the block diagram 
are elements with names INI and IN2, so that the total number of elements 
in the block diagram is six. All of the information concerning these 
elements is contained within the ELEMENTS statement shown in Figure 3. 

When elements are present, they must be connected together to form the 
block diagram. This is done using a CONNECT statement. One important 
rule that must not be violated xjhen making connections is that redun- 
dant connections are to be avoided. The same path must not be re- 
traced. Since there is great flexibility in constructing CONNECT state- 
ments, the user should always strive to keep it simple for "bookkeeping" 
purposes. The first CONNECT statement in Figure 3 connects the output 
of adder ADI to an output terminal AD10UT. Note that the dollar sign ($) 
terminates a string of connections. The next CONNECT statement (denoted 
by the colon) traces the feed forward path and the lower feedback path, 
all in one statement. This is done strictly for the sake of convenience 
and speed of writing. This connection could have been effected using 
many separate CONNECT statements, if desired. The next CONNECT statement 
connects the output of integrator INI with output terminal XDOT. The 
next statement traces the upper feedback path and the last statement con- 
nects integrator IN2 with output terminal X. CONNECT statements are easy 
to construct and should present no difficulty for the user. More informa- 
tion on CONNECT statements can be found in Section II of this Manual. 
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When the problem i£ simulated in EQUATION form, ELEMENTS and CONNECT 
statements are unnecessary. All that is required is the EQUATION oper- 
ator as shown in Figure 4. Each equation is terminated by a ($) sign. 

The colon repeats the EQUATION operator for each equation. 

The END$ statement terminates the Description Module. 

Simulation Module 

The SIMULATE statement name MASS SPRING DAMPER references the model 
being simulated. The name is the same as the main MODEL name given in 
the Description Module. The excitation function, sin t, is impressed 
upon the input terminal U using the EXCITE statement whose format is 
given in the Table of Standard Excitation Functions. The integration 
method selected to solve this problem is Runga-Kutta with a step size of 
0.01. This information is recorded in the INTEGRATE statement. The 
STOP IF statement terminates the simulation after five time units. The 
initial condition X(0) = 20.0 is specified in the INITIALIZE statement. 

In Figure 3, the initial condition is placed on the integrator whereas 
in Figure 4, X is initialized directly. 

The END$ statement terminates the Simulation Module. 

Post Processing Module 

Since we wish to examine the system input and outputs, the PRINT 
statement contains all of the names listed in the INPUTS and OUTPUTS 
statement which are to be printed put. TIME will automatically be listed 
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in the first column of each page of printout. The printing will occur 
at each integration STEP as indicated in the SAMPLE statement. If the 
SAMPLE statement is omitted, printing will occur at each integration 
step. 

The END$ statement terminates the Post Processing Module. 
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G. MARSYAS MATH REFERENCE 

When using the EQUATION option in the Description Module, mathe- 
matical expressions may arise which allude to a quantity called a MATH 
REFERENCE. The MARSYAS MATH REFERENCE library consists of the Standard 
Elements, user-defined elements (DEVICES) and Standard Excitation Func- 
tions. The form which these MATH REFERENCES take may be quite complex, 
since they may consist of mathematical expressions, mathematical opera- 
tions or other MATH REFERENCES. 

Suppose an equation contains a forcing function f(x,y,t) which can 
be represented by a Standard Excitation Function. The general form of 
the Standard Excitation Functions, as shown in the Table of Standard 
Excitation Functions at the back of the Manual is as follows: 

excitation function mnemonic (parameters) $ 

For use in an EQUATION as a MATH REFERENCE, this form must be 
altered slightly to include TIME or an equivalent expression for time 
as 1 

excitation function mnemonic (parameters. A) $ 

where A is either the word TIME or a mathematical expression from which 
time can be computed. Usually, however, the user will simply insert the 
word TIME in place of A. (See EXAMPLE (2) in EQUATION, Part II A.) The 
user then associates each of the parameters required by the Standard 
Excitation Function with the variables in f(x,y,t) as shown in EXAMPLE (1). 
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EXAMPLE (1) 

Suppose an equation contains a forcing function 5X sin(2X^ + 3Yt). 
Using the Table of Standard Excitation Functions at the back of the 
Manual, the MATH REFERENCE is coded as follows: 

FSIN(5. * X, 3. * Y, 2.0 * X ** 2, TIME) 

where d x = 5. * X, d 2 = 3. * Y, and d 3 =2. * X ** 2. 

Next, suppose that an EQUATION contains a term which is to be repre- 
sented by a Standard Element. This situation might arise when using 
EQUATION in conjunction with a block diagram. In this case, the MATH 
REFERENCE is represented in much the same way as with the Standard 
Excitation Function except that we now use Standard Elements in place 
of Standard Excitation Functions. The general form of the Standard 
Elements as shown in the Table of Standard Elements at the back of the 
Manual is as follows: 

element mnemonic (parameter s) $ 

For use as a MATH REFERENCE, this form must be altered slightly 
to include a mathematical expression(s) for each input of the element, 
as well as an integer indicating to which element output terminal the 
EQUATION applies, as 

element mnemonic (parameters, M^, M 2 , ..., N) 

where M^ is a mathematical expression for each element input terminal 
which may involve derivatives, other mathematical expressions or other 


31 



/ 


I G. MARSYAS MATH REFERENCE 

MATH REFERENCES. N is the integer indicating to which output terminal 
the EQUATION applies. 

If the element has only one output, the integer N may be omitted. 
The user assigns numerical values to the indicated parameters (if any) 
as required in the usual use of Standard Elements. See EXAMPLE (2). 

Note in EXAMPLE (2) that, since there are no parameters associated 
with an integrator element and only one output, both the parameters and 
the integer N are omitted from the parenthesis. Additionally, only one 
mathematical expression is needed since the integrator has only one in- 
put. Notice that either method of coding 0UTPUT1 is legitimate since 
both are equivalent mathematically. 

EXAMPLE (2) 

Suppose a block diagram contains two integrators and one constant 
multiplier as shown in the sketch below. Using the mnemonic for an 
integrator, IN, as given in the Table of Standard Elements at the back 
of the Manual, an equation for 0UTPUT1 using a MATH REFERENCE is coded 
as shown below. 
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EQUATION: Z" = U - 3.2 * Z' $ 

: OUTPUTl = IN(Z") $ 

or 

EQUATION: OUTPUTl = IN(U - 3.2 * Z *) $ 
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II. MODULAR STRUCTURE OF MARSYAS 


A. Description Module 

The function of the Description Module is to provide the MARSYAS 
processor with a description of the structure of the system being simu- 
lated. The system might be a mathematical representation of a complete 
Space Shuttle, its engines or any one of its subsystems. The mathematical 
representation could be in the form of a block diagram, a series of block 
diagrams, a set of differential equations or a mixture of block diagrams 
and equations. 

A Description Module is a MARSYAS representation of a complete or 
partially complete mathematical system. A MARSYAS program may contain 
several Description Modules. A model is completely described in a 
Description Module. A submodel is a model wholly contained within another 
model. If a model contains one or more submodels, then each submodel's 
structure must be described in a separate Description Module. 

When a model is represented in block diagram form, its structure is 
specified by its INPUTS, OUTPUTS, ELEMENTS, SUBMODELS and connections. 

An input is any point in a model where a driving function is applied and 
an output is any point where the behavior of the model is to be examined. 
An element is any device whose input (s) and output (s) can be related by 
algebraic or differential equations or based on logic. Each element 
type, such as an adder or integrator, has a mnemonic, viz., AD and IN, 
respectively. A library of Standard Elements is part of the MARSYAS 
system and is given in the Table of Standard Elements at the back of 
the Manual. A FORTRAN subprogram, called a DEVICE, can be constructed 
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* 


by the user to simulate any element not found in the Table of Standard 
Elements. Some elements contain parameters which can be numerical con- 
stants or named constants. If the parameter constants are named, their 
values are specified with PARAMETERIZE statements. The interconnections 
among the system inputs, elements, submodels and outputs are specified 
with CONNECT statements. Connect statements can be modified without 
altering the original coding by using DISCONNECT statements. Alternate 
names for models, elements, inputs, outputs or parameters can be assigned 
using the NAME statement. There are no reserved names (words) within the 
MARSYAS system. 

When a model is represented by a collection of differential equations 
the statement operators ELEMENTS, PARAMETERIZE, SUBMODEL, CONNECT, 
DISCONNECT and NAME are not used. The EQUATION operator is used when 
specifying equations. Equations may contain mathematical expressions 
which consist of MARSYAS names, mathematical operators, numerical con- 
stants and MATH REFERENCES. MATH REFERENCES are discussed in Section I G 

The order in which the above statements appear within a model 
description is immaterial except that MODEL and END must be the first 
and last statements, respectively. A detailed description of all of the 
statement operators used in the Description Module appears in alphabeti- 
cal order on the following pages. Figure 5 is a list of these operators 
with a brief description of their function. 
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DESCRIPTION MODULE 


OPERATOR 

CONNECT 

DEVICE 

DISCONNECT 

ELEMENTS 

END 

EQUATION 

INPUTS 

MODEL 

NAME 

OUTPUTS 

PARAMETERIZE 

\ 

i 

SUBMODEL 


FIGURE 5. STATEMENT 


FUNCTION 

Specifies the interconnection of ele- 
ments, input and output terminals and 
submodels within the model. 

Defines a non-standard element. 

Cancels connections specified by a 
previous CONNECT statement. 

Specifies the mnemonic, names and 
parameters of elements in a model. 

Marks the end of a Description Module. 

Specifies a model in terms of 
equations. 

Assigns names to the model input 
terminals . 

Marks the appearance of a Description 
Module and names the model. 

Assigns alternate names to units already 
defined. 

Assigns names to the model output terminals 

Assigns a numerical value to a parameter 
which was given a name instead of a num- 
ber in an ELEMENTS statement. Also used 
to change the value of a parameter as- 
signed a numerical value by an ELEMENTS 
or PARAMETERIZE statement. 

Assigns a name to a submodel and its 
associated inputs and outputs. 

Allows comments for the user. 


OPERATORS USED IN THE DESCRIPTION MODULE 
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/ 


COMMENT 


FORMAT - 


* message or blanks $ 

COMMENTARY - 

The comment statement may appear anywhere in any module and may 
extend over as many source cards as are necessary. The source cards 
are printed exactly as read without the editing common to all other 
MARSYAS statements. The asterisk can appear in any card column. 

EXAMPLE 

* This is a comment card $ 


i 

i 
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CONNECT 

FORMAT - 

CONNECT: namel, name2, name 3, etc. $ 

COMMENTARY - 

The CONNECT statement is used to describe the signal paths among 
the elements, submodels and input and output terminals of the system. 
The direction of the signal flow is implied by the order in which the 
arguments of the connection statement are given. CONNECT statements 
must always begin as follows: i 

1. at an element output 
or 

2. at a system input terminal 
CONNECT statements must always end as follows: 

1. at an element input , 
or 

2. at a system output terminal 

Since a system input terminal is used to apply a signal, it must 
be the first to occur in the CONNECT statement of which it is a part. 
Similarly, since a system output name can only receive a signal, it can 
only appear as the last name in a string of CONNECT statements. Redun- 
dant connection paths are not allowed. 
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CONNECT (Cont'd.) 

If a signal path goes through a submodel, it is only necessary to 
specify the connections to and from the submodel. It is not necessary 
to repeat the connections within the submodel which have previously been 
described in that submodel's Description Module. 

Every signal path which occurs within the model must be described 
by a CONNECT statement. If a signal path branches, then the new path 
is started at the point of branching. The CONNECT statement is not used 
with the EQUATION statement. 


EXAMPLE 


INPUT1 


oi INI 


ADI 


0UTPUT2 


The CONNECT statement for the above simple circuit is as follows: 
CONNECT: INPUT1, INI, ADI, 0UTPUT2$ 


EXAMPLE 

Consider the following simple circuit: 
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CONNECT (Cont'd.) 

The CONNECT statement reads as follows: 

CONNECT: INPUT5, A, B, C, D, OUTPUT? $ 

: B, E, - A$ 

There are several permissible variations of this CONNECT statement 
such as 

CONNECT: INPUTS, A, B, E, -A 

: B, C, D, OUTPUT? $ 

The intermediate dollar sign(s) may be omitted when using the colon 
to repeat a statement operator as above. 

Note that the CONNECT statement allows the use of a negative sign 
on element inputs only. For example, -INPUT5 would not be permitted in 
the above CONNECT statements. 

When elements with multiple inputs and outputs are connected, the 
element terminal numbers must be carefully noted when writing the CONNECT 
statement. 
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II A. DESCRIPTION MODULE 

CONNECT (Cont'd.) 

EXAMPLE 

Consider the following connections of system INPUT5 , A1 and B1 
with ELEMENTS, ML1 and IN4 and the connection of the output of element 
ML1 with the system output, Cl. 



The CONNECT statement could be coded as follows: 

CONNECT: Bl, IN4, j l#MLl, ML1#1 j , Cl$ 

: Al, 2#ML1$ 

The two terms in the dotted box may be combined into one for ease 
of writing. Thus l#MLl#l is equivalent to writing 1#ML1, ML1#1. 

The above connect statement reads as follows: "Connect input Bl 

to element IN4 and element IN4 to input #1 of element ML1 and connect 
output #1 of element ML1 to output Cl. Connect input Al to input #2 of 
element ML1." Note that the second CONNECT statement, CONNECT: Al, 2#ML1$, 

does not connect past input #2 of ML1, thus avoiding a redundant connection. 
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CONNECT (Cont'd.) 

Should the user desire, he can give names to the terminals of an 
element using a NAME statement, then CONNECT using element names instead 
of terminal numbers (see NAME) . 

EXAMPLE 

Suppose that we wish to assign the names INPUT1, INPUT2 , and OUTPUT 
to the three terminals of ML1 in the previous example. The coding for 
the NAME and CONNECT statements is as follows: 

NAME: 1#ML1 , INPUT 1$ : 2#ML1, INPUT2$ : MLl#l, OUTPUT$ 

CONNECT: Bl, IN4, INPUT1, OUTPUT, Cl$ 

: Al, INPUT2$ 

When elements with multiple inputs and a single output or a single 
input and multiple outputs are connected, additional simplification of 
the coding is possible as shown in the EXAMPLE below. 
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CONNECT (Cont'cQ 


EXAMPLE 


The three elements with names A, B and C are connected as shown 
in the sketch below. 


INI 

IN2 

IN3 



The CONNECT statement can be coded as follows: 

CONNECT: INI, 1#A#1, B, 1#C#1, 0UT1 $ 

: IN2, 2#A $ : IN3, 3#A $ 

: C#2, 0UT2 $ : C#3, 0UT3 $ 

The first CONNECT statement can also be written as follows: 
CONNECT: INI, 1#A, B, C#l, OUTl $ 

This shortened form of the CONNECT statement is possible since A 
has a single output, B has a single input and output and C has a single 
input . 
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DEVICE 


FORMAT - 

DEVICE: a FORTRAN name, number of element input terminals, 

number of element output terminals, number of element 
parameters $ 


COMMENTARY - 

The DEVICE statement is used to create a new Element, coded in 
FORTRAN, which is treated as though it were a Standard Element. This 
is done by specifying the mnemonic, as identified by the FORTRAN name, 
and referenced in a MARSYAS ELEMENTS statement, and the number of inputs, 
outputs and parameters the new element is to have. A FORTRAN subroutine 
must immediately follow the DEVICE statement and is used to add to the 
MARSYAS library a temporary program capable of calculating the output of 
the new element from its input values. A DIMENSION statement is required 
in the subroutine even if there are no arrays needed within the subrou- 
tine. 

The subroutine name referenced in the FORTRAN SUBROUTINE statement 
is the same FORTRAN name referenced in the DEVICE statement. This is the 
element mnemonic assigned in the appropriate Description Module ELEMENTS 
statement. A maximum of four arguments in parenthesis may folloxtf the sub- 
routine name. These arguments are arbitrary real FORTRAN array names for 
the element inputs, outputs, parameters and time, in that order. If time 
is not needed in the subprogram, it may be omitted from the argument list. 
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DEVICE (Cont’d.) 

If time is used, however, its value must not be altered or otherwise 
manipulated within the subprogram. Since the element inputs, outputs 
and parameters are treated as arrays within the FORTRAN subprogram, a 
DIMENSION statement dimensioning the arrays must appear within the sub- 
routine. Real variable array names must be given to the dimensioned 
arguments. The element name referenced in the MARSYAS ELEMENTS state- 
ment is a MARSYAS name and, hence, is not restricted to six characters 
as are FORTRAN names. 


EXAMPLE (1) 

The user decides to construct an element with MARSYAS name ALPHAONE 
for use in Description Module model DELTA. The element, having a FORTRAN 
mnemonic name ALPHA, has three inputs, two outputs and eight parameters. 
The arbitrary real array names for the element inputs, outputs, parameters 
and time are A, B, C, and D, respectively. The coding is as follows: 


CO 

< 

60 


G 

CO 

*H 

1 

O 

a 


DEVICE: ALPHA, 3, 2, 8 $ 


SUBROUTINE ALPHA (A, B, C, D) 
DIMENSION A(3) , B(2), C(8) 


• 

• • ■ 

• 

• 

RETURN 

• • 

• 

END 



MODEL: 

• 

DELTA 

• • 

• 

• 

• • 

• 


ELEMENTS: ALPHA, ALPHAONE (1.1, -3.0, 0.0, 6.2, 

8.0, 2.1, 3.4, -7.6) $ • 

• • • # 

END$ 
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DEVICE (Cont'd.) 


EXAMPLE (2) 

Shown below is the FORTRAN coding for a relay device which i 
available in the Standard Elements library of MARSYAS. Following it 
the MARSYAS coding referencing this device. 



DEVICE: RELAY, 1, 1, 2$ 

i SUBROUTINE RELAY (X, Y, Z) 

DIMENSION X(l), Y(l) , Z(2) 

YBREAK = Z(l) 

<j m XBREAK = Z(2) 

g £ IF (ABS X(l) . LT. XBREAK) Y(l) = 0.0 

go IF (X(l) . GE, XBREAK) Y(l) = YBREAK 

IF (X(l) .LE. -XBREAK) Y(l) = -YBREAK 

RETURN 

END 

• • • • • 

MODEL: TIM$ 

INPUTS: ... $ 

OUTPUTS: ... $ 

ELEMENTS: RELAY, RELAY1 (YBREAK, XBREAK) $ 

CONNECT: ... $ 

PARAMETERIZE: RELAY1, (YBREAK, 2.0, XBREAK, 1.0)$ 

END$ 


Note that the FORTRAN name for the DEVICE and for the input , 
and parameter array names must be real variable names. 


s 

is 


output 
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DEVICE (Cont'd.) 


Should the need arise, the DEVICE statement can be constructed in 
a general manner by substituting the letter N for any or all of the num- 
ber of element inputs, outputs and parameters. Then the value (s) of N 
is specified in the ELEMENTS statement as shown in EXAMPLE (3) below, 
following by the appropriate element parameters. 


EXAMPLE (3) 

If in EXAMPLE (2) the number of element inputs, outputs and param- 
eters are to be varied, the coding would appear as follows: 


DEVICE: ALPHA, N, N, N $ 

SUBROUTINE ALPHA (A, B, C, D) 
DIMENSION A(3) , B(2), C(8) 


RETURN 

END 


MODEL: 


MODEL: 


PSY$ 


BETA$ 


ELEMENTS: ALPHA, ALPHA1(2, 1, 4, '3.2, 1.1, 6.4, 9.3)$ 

END$ 


ELEMENTS: ALPHA, ALPHA15(1, 2, 3, 0.0, 2.4, 6.1)$ 

END$ 
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DEVICE (Cont'd.) 

Note in EXAMPLE (3) that the first three arguments in the element 
parameter string identify the values of N to be inserted in the device 
statement. The first value of N being the number of element input 
terminals, the second value of N being the number of element output 
terminals and the third value of N being the number of element parameters. 
The remaining numbers in parenthesis are the actual values of the ele- 
ment parameters. 

It is permissible to place DEVICE statements, and associated sub~ 
routines anywhere within the MARSYAS program after the BEGIN statement 
and before the Simulation Module provided they appear prior to their 
first reference in an ELEMENTS statement. It is recommended, however, 
that all DEVICE statements be placed at the very beginning of the MARSYAS 
deck, after the BEGIN statement. This serves the purpose of isolating 
the source FORTRAN code from the MARSYAS language statements. 
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DISCONNECT 


FORMAT - 


DISCONNECT: namel, name2, etc. $ 

COMMENTARY - 

The DISCONNECT statement is used to cancel already defined signal 
paths among the elements, submodels and input and output terminals. The 
statement is written in exactly the same manner as for the CONNECT state 
ment. The DISCONNECT must appear after a CONNECT statement which estab- 
lished the connection which DISCONNECT is now to break. The DISCONNECT 
statement is not used with EQUATION. 

EXAMPLE 


If a CONNECT statement is given as 
CONNECT: A, B, C, D $ 

a DISCONNECT statement might be 

DISCONNECT: B, C $ 

The above two statements are now equivalent to the following statements: 


CONNECT: A, B $ 

CONNECT: C, D $ 
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ELEMENTS 


FORMAT - 

ELEMENTS: mnemonic, name (parameters, if required) $ 

COMMENTARY - 

MARSYAS provides the user with a larger number of operational 
devices called ELEMENTS, similar to those available on an analog computer, 
but far more numerous and versatile. These elements include items such 
as adders, integrators, multipliers, transfer functions, etc. A Table of 
Standard MARSYAS ELEMENTS appears at the back of the Manual. The Table 
lists the elements according to type or class, gives the element block 
diagram symbol, the number of inputs and outputs, the element mnemonic, 
the input-output relation and the appropriate parameters, if any. The 
user is not restricted to the "standard" set. Non-standard elements may 
be defined by the user through the use of the DEVICE statement. These 
new elements are treated in exactly the same manner as the original set 
furnished by the system. The ELEMENTS statement is not used when the 
EQUATION statement is used. 

EXAMPLE 

ELEMENTS: AD, ADDER1 , ADDER2 , ADDER3 $ 

ELEMENTS: CM, CMl(l.O), CM2 (-3.6) $ 
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ELEMENTS (Cont'd.) 

The constant multiplier (CM) elements above specify the parameter 
constants 1.0 and -3.6. If the user desires, dummy variables may be sub- 
stituted for these values and later specified in PARAMETERIZE statements. 

An example of using parameter constants in an ELEMENTS statement 
is as follows: 

ELEMENTS: CM, CPU (A) , CM2 (B) $ 


PARAMETERIZE: CM (A, 1.0) $ 

: CM2 (B, -3.6) $ 


See CONNECT for additional comments on elements. 
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END 


FORMAT - 


END$ 


COMMENTARY - 

The END statement is used to close a DESCRIPTION MODULE. The 
END statement must be the last statement in a DESCRIPTION MODULE. 


EXAMPLE 

END$ 
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EQUATION 

FORMAT - 

EQUATION: Differential or algebraic equation $ 

COMMENTARY - 

The EQUATION statement operator is designed to provide the user 
with great flexibility in specifying his model. The model can be 
described in terms of a system of differential and algebraic equations 
with almost no restrictions on their structure. The one main restriction 
is that a "solution scheme" of the set of equations must exist. This is 
discussed in detail in Section IV C. of the Manual. 

The ordinary differential equations which MARSYAS solves may be 
linear or nonlinear and of any degree. Since the independent variable 
is always time in the MARSYAS system, all differentiation is assumed to 
be with respect to TIME. Direct differentiation of the variable name 
TIME or of system INPUT or OUTPUT names is not permitted. The output vari- 
bles may be differentiated .by assigning to them alternate names. An 

apostrophe is used to indicate differentiation with respect to time. 

9 9 

Thus d^x/dt^ would be represented as X". 

There are no restrictions on the form of the mathematical expressions 
contained within the equations, except that the expression itself must not 
be differentiated. Expressions may employ any of the standard arithmetic 
operations of addition (+) , subtraction (or minus sign) (-), multiplica- 
tion (*) , division (/) and exponentiation (**) . Balanced parentheses may 
be used where needed. See EXAMPLE (1). 


r 
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EQUATION (Cont'd.) 

Mathematical expressions may consist of MARSYAS names, mathematical 
operators, numerical constants and MATH REFERENCES. The use of the 
MATH REFERENCE adds an additional dimension to the use of the EQUATION 
option in that the user is permitted to take shortcuts when constructing 
equations, especially when the equations are used in conjunction with 
block diagrams. See MATH REFERENCE, Section I G, 

The block diagram and equation modes may be mixed when the system 
main MODEL, described in block diagram form, contains SUBMODELS, some or 
all of which contain equations. The structure of each submodel is, of 
course, described as a separate MODEL. 

When a MODEL contains an EQUATION statement, the statement operators 
ELEMENTS, PARAMETERIZE, SUBMODEL, CONNECT, DISCONNECT and NAME are not 
necessary and, therefore, not used. The MODEL statement must contain the 
word EQUATION to alert the MARSYAS processor that equations are used in 
the MODEL (see MODEL). 

Many times, the EQUATION statement is used, in conjunction with the 
INPUTS and OUTPUTS statements, to build a model of a system. When used 
in this manner, the EQUATION operator presents a mathematical relation 
which determines the value of the system outputs from the value of the 
inputs and TIME. These relations may involve an arbitrary number of 
intermediate variables. 

An example of using MARSYAS in this manner is illustrated in Fig- 
ure 4, Section I F. Note that each equation in Figure 4 is preceded by 


/• 


55 



J 


II A. DESCRIPTION MODULE 

O 


EQUATION (Cont'd.) 


the EQUATION statement (or its colon equivalent) and followed by the 
dollar sign ($) terminator. 

The number of INPUTS and OUTPUTS names for a MODEL is arbitrary. 
There need not be any INPUTS, but there must be at least one output 
terminal when using EQUATION. Thus, the INPUTS statement may be omitted 
but the OUTPUTS statement must appear. 

Initial conditions on the equations are coded in the Simulation 
Module as noted in the INITIALIZE statement in Section II C. of this 
Manual . 

EXAMPLE (1) 


Code the following homogeneous differential equation using the 
EQUATION operator. 



dx 

dt 


+ 


cx 


0.0 


EQUATION: X" + X' + C*X = 0.0$ 
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EQUATION (Cont'd.) 


EXAMPLE (2) 


Code the following non-homogeneous differential equation using the 
EQUATION operator and the MATH REFERENCE. 



+ 


dx 

dt 


+ cx 


5 sin 2 t 


EQUATION: X" + X 1 + C*X = FSIN (5.0, 2.0, 0.0, TIME) $ 


An example of using the EQUATION operator with a MARSYAS transfer 
function element converted to a MATH REFERENCE is shown in EXAMPLE (3) 
of INITIALIZE in Section II C. 
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INPUTS 


FORMAT - 


INPUTS: namel, name2, etc. $ 

COMMENTARY - 

The INPUTS statement is used to give identifying names to the 
points in the model (input terminals) at which driving functions are 
to be introduced. The INPUTS statement is optional if driving func- 
tions are not needed in the model. 

EXAMPLE 


INPUTS: INI , IN2, IN3 $ 

or 

INPUTS: INI $ 

: IN2 , IN3 $ 
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MODEL 


FORMAT - 


MODEL: name$ or MODEL : name, EQUATION $ 


COMMENTARY - 

The MODEL statement is used to give an identifying name to the 
model'’ or submodel being described. Within a Description Module, all 
model names must be unique. MODEL must be the first statement of a 
Description Module and END the last. The name referenced is the user- 
specified model name. 

EXAMPLE 

MODEL: SATURN5 SIMULATOR$ 

MODEL: X-LOOP GYRO$ 

MODEL: ALPHA$ 

If the EQUATION option is chosen, then the MODEL statement in- 
cludes the word EQUATION as shown in the example below. 

EXAMPLE 

MODEL: DELTA, EQUATION$ 
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NAME 


FORMAT - 


NAME: original name, new name $ 

COMMENTARY - 

The NAME statement is used to assign an additional MARSYAS 
name to an element, parameter, system input terminal, or system output 
terminal which is defined elsewhere through an ELEMENTS, PARAMETERIZE, 
INPUTS, or OUTPUTS statement, respectively. In addition, it may be used 
to assign a new name to a particular terminal of an element or its asso- 
ciated parameters. The NAME statement is not used in connection with the 
EQUATION statement. 

EXAMPLE 


NAME: BILL, WILLIAM$ 

This NAME statement assigns the new name WILLIAM to BILL. The 
name BILL is not erased from the MARSYAS program and the names WILLIAM 
and BILL may be freely interchanged. 

The NAME statement can also be used to assign names to the 
input(s) and output(s) of elements. See CONNECT for an example. 

It is illegal to attempt to change a name previously defined by 
a NAME statement with a new NAME statement. 
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OUTPUTS 


FORMAT - 


OUTPUTS: namel, name2, etc. $ 

COMMENTARY - 

The OUTPUTS statement is used to give identifying names to the 
points in the model (output terminals) at which the system is to be 
examined. There must be at least one output terminal in the main model. 

EXAMPLE 


OUTPUTS: OUT1 , 0UT2, 0UT3$ 

or 

OUTPUTS : 0UT1$ 


0UT2 , 0UT3$ 
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PARAMETERIZE 


FORMAT - 


or 


PARAMETERIZE: element name (parameter number or name, numerical 

value) $ 

PARAMETERIZE: parameter name (numerical value) $ 


COMMENTARY - 

The PARAMETERIZE statement is used to give numerical values to 
specific MARSYAS parameters or to change parameter values previously 
assigned. Parameters are identified through the names given in ELEMENTS 
and NAME statements. The PARAMETERIZE statement is not used with 
EQUATION. 


EXAMPLE 

Suppose the following is referenced in an ELEMENTS statement 

ELEMENTS: TF, TF1(2, 2.0, 1.0, PHY, 1.8, 8.5, BETA) $ 

Then the PARAMETERIZE statement might read 

PARAMETERIZE: TF1 (PHY, 6.0)$ : TF1 (BETA, -1.0E-2)$ 

or 

PARAMETERIZE: PHY (6.0)$ : BETA (-1.0E-2)$ 
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PARAMETERIZE (Cont'd.) 

The PARAMETERIZE statement can be used to alter the value of a 
parameter previously specified in an ELEMENTS statement. The user may 
wish to use this method to effect a parameter change rather than recode 
the original ELEMENTS statement, since it requires less work. If an ele- 
ment contains several parameters, the PARAMETERIZE statement effecting a 
change in parameters can be shortened by specifying the position of the 
parameter within the parentheses, along with its revised value. See 
CHANGE in Section II C. 

EXAMPLE 


For the ELEMENTS statement shown below, change the value of the 
fifth parameter from 1.8 to 15.2. 

ELEMENTS: TF, TF6(2, 1.3, 2.1, 7.4, 1.8, 3.4, 9.7) $ 

PARAMETERIZE: TF6(5, 15.2)$ 


The ELEMENTS then appears to the system as if it had been origi- 
nally coded as 

ELEMENTS: TF, TF6(2, 1.3, 2.1, 7.4, 15.2, 3.4, 9.7) $ 


63 



II A. DESCRIPTION MODULE 


j 


SUBMODEL 


FORMAT - 


SUBMODEL: name 


INPUTS: a, b, c, d, etc; OUTPUTS: h, i, j, k, etc$ 


COMMENTARY - 


A submodel is a model wholly contained within another model. The 
SUBMODEL statement is needed to give identifying names to the submodel and 
its input and output terminals. These terminal names will be used to 
CONNECT the submodel to the main model. Submodels can be imbedded within 
submodels. The structure of each submodel must be described in a Description 
Module separate from the one in which it is imbedded. 


EXAMPLE 


MODEL: GAMMAS 



64 





II, 'A. DESCRIPTION MODULE 


SUB MODEL (Cont'd.) 


The model GAMMA contains one submodel and an integrator. It 
could be cocied as follows: 


o«o*o 


i 

r~ 


MODEL: 

GAMMA $ 




INPUTS: 

INPUT1 , INPUT2$ 

CD 

rH 



OUTPUTS: 

OUTPUT 1 , OUTPUT 2, 0UTPUT3$ 






£ 

<g 


ELEMENTS: 

INTEGRATOR? 

c 

o 

•H 

■U 

o 

LJ 

SUBMODEL: 

PHY; INPUTS: Al; OUTPUTS: Bl, B2? ] 

P< 

•H 

U 

o 

o 

IH 


CONNECT: 

INPUT 1, Al$ 

w 

<D 

P 




Bl, OUTPUT 1$ 





B2, OUTPUT 2$ : INPUT2, INTEGRATOR, 0UTPUT3 $ 

'l 

r 


END? 






• • • 

/ 

k 

1 — 

— ““ — 

PHY $ "j 

O 

r 

L_ 

MODEL: 

i — 1 

3 

jg 



INPUTS: 

Al$ 

£ 

CM 


OUTPUTS : 

Bl, B2$ 

o 




•H W 
4J O 

•H 


ELEMENTS : 

n 

> Describes the configurations of PHY 

O 

W 



CONNECT: 

* ) 
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SUBMODEL (Cont'd.) 

NOTE: In describing the submodel, the input and output terminal 

names of the submodel need not be the same as those defined in the main 
model. However, there must be a one-to-one correspondence between the 
number and order of the terminals specified in both models' INPUTS and 
OUTPUTS statements. 

MARSYAS names used within a particular MODEL must be unique. 
However, names which are used in a submodel can be repeated in the main 
MODEL. 
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II B. MODIFICATION MODULE 
B. Modification Module 

The Modification Module is used to alter an already existing model 
by adding or deleting elements, input or output terminals, and submodels 
and/or by changing the interconnections among these items. A MARSYAS 
program may contain several Modification Modules. 

Figure 6 is a list of statements which can appear in the Modification 
Module with brief comments as to their function. A detailed explanation 
of a few of those statements, not previously discussed, is given in the 
following pages. Explanations of the remaining statements appear in the 
Description Module portion of this Manual. The order in which the state- 
ments appear is immaterial except that MODIFY must be first and END last. 
DISCONNECT, DELETE and SUBSTITUTE can be used to alter the Description 
Module statements which defined the original connection (CONNECT) which is 
now to be broken or the original element (ELEMENTS) which is now to be 
deleted or substituted. DEVICE may appear at any point in the program pro- 
vided it is used before the first appearance of the new element type in an 
ELEMENTS statement. The INPUTS and OUTPUTS statements are used to define 
additional input and output terminals, respectively. The SUBMODEL, ELEMENTS 
and DEVICE statements are used to add additional submodels, elements and 
non-standard elements, respectively. All additional terminals, if any, must 
be connected within the Modification Module. The PARAMETERIZE and NAME 
statements have the same functions, as described within the Description 
Module. When the EQUATION option is used in the Description Module, there 
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can be no Modification Module since no ELEMENT'S, element names, CONNECTS 
or parameters will exist there. 

If the system contains several Description and Modification Modules, 
the only requirement for ordering the Modules within a MARSYAS program 
is that the Modification Module must follow the Description Module it 
modifies. For example, a MARSYAS program might contain a group of three 
Description Modules followed by a group of three Modification Modules. 
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OPERATOR 

CONNECT 

DELETE 

DEVICE 

DISCONNECT 

ELEMENTS 

END 

INPUTS 

MODIFY 

NAME 

OUTPUTS 

PARAMETERIZE 

SUBMODEL 

SUBSTITUTE 

* 

FIGURE 6. 


FUNCTION 


Same functions as in Description Module. 

Deletes an already specified element or 
input or output terminal, and removes 
all connections to or from that element 
or terminal. 

Defines a non-standard element. 

Same functions as in Description Module. 

Adds additional elements to the model. 

Marks the end of a Modification Module. 

Adds additional input terminals to the 
model . 

Marks the appearance of a modification 
module and can give a new name to the 
modified model. 

Same functions as in Description Module. 

Adds additional output terminals to the 
model . 

Assigns or changes numerical values of 
parameters of elements present in 
either a Description or Modification 
model . 

Adds additional submodels to the model. 

Substitutes one element for a previously 
defined element. The number of inputs 
and outputs must be the same for both 
elements . 

Same functions as in Description Module. 


STATEMENT OPERATORS USED IN THE MODIFICATION MODULE 
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DELETE 


FORMAT - 


DELETE: namel, name2 , name3 , etc,$ 

COMMENTARY - 

The DELETE statement removes an element, input terminal or 
output terminal from the model being modified. All connections asso- 
ciated with the deleted items are broken. The user must reconnect 
(using CONNECT statements) , the elements and terminals which are 
affected by the DELETE statement. 

EXAMPLE 


If the unmodified model ELEMENTS statement is 
ELEMENTS: AD, ADI , AD2 , AD3$ 

and a Modification Module DELETE statement is 
DELETE: AD2$ 

then the effect is that the original model would look as if it had 
been coded as 

ELEMENTS: AD, ADI , AD3$ 
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DELETE (Cont'd.) 


The user must reconnect the deleted connections associated with 
the deleted element, AD2. 


\ 

! 
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END 


FORMAT - 


END$ 


COMMENTARY - 

The END statement is used to close a Modification Module. 
END statement must be the last one in a Modification Module. 


EXAMPLE 


END$ 


I 
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FORMAT - 


MODIFY 


MODIFY: model name, new model name (optional) $ 


COMMENTARY - 

The MODIFY statement is used to identify the model which is to be 
modified and optionally to give the modified model a new name. If a 
new name is not given, the modified model has the same name as the ori- 
ginal model. MODIFY must be the first statement of the Modification 
Module. 


EXAMPLE 

MODIFY: ALPHA$ 

MODIFY: ALPHA, BETA$ 

The first statement identifies the model which is to be modified. 
The second statement may be used in place of the first if the user de- 
sires to retain the original model, ALPHA, and create a new model BETA 
which will be a modified version of ALPHA. 
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SUBSTITUTE 


FORMAT - 


SUBSTITUTE: original element name , replacement element name $ 

COMMENTARY - 

The SUBSTITUTE statement is used to effect the complete replace- 
ment of one element by another element. The elements do not have to be 
of the same type as long as each element has the same number of input 
and output terminals. For example, a power function could be substituted 
for an integrator* since both have one input and one output. All connec- 
tions to and from the original element are replaced by the connections 
to and from the substitute element. 

EXAMPLE 

If the original model contained the following statements, 

ELEMENTS: CM, CM1$ 

: AD, ADI, AD2$ 

CONNECT : ADI, CML, AD2$ 

and the Modification Module contained the statements, 

SUBSTITUTE: CML, INL$ 

ELEMENTS: IN, IN1$ 
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SUBSTITUTE (Cont'd.) 

then the effect is that the original model would look as if it had been 
coded as 


ELEMENTS: AD S ADI, AD2$ : IN, IN1$ 

CONNECT: ADI, INI, AD2$ 


Note: It is up to the user to determine whether a particular sub- 


stitution is physically meaningful. 



II C„ SIMULATION MODULE 


C. Simulation Module 

The Simulation Module completes the description of the system whose 
structure is defined in a Description Module (s). A MARSYAS program can 
contain only one Simulation Module. In this Module, the user specifies 
the INPUT excitation function (s), the numerical integration scheme to be 
used in the problem solution and any required initial conditions. The 
conditions for stopping the simulation are also specified in the Simula- 
tion Module. 

A sequential series of simulations can be run without resubmitting 
the MARSYAS deck by using CHANGE statements. These permit the user to 
alter one or more of the following statements without altering the origi- 
nal coding: EXCITE, INITIALIZE, STOP IF, TERMINATE IF and PARAMETERS . 

Excitation functions are impressed upon the model input terminals 
using EXCITE statements. A Table of Standard Excitation functions is 
part of the MARSYAS system and is shown at the back of this Manual. 
Additional excitation functions may be defined by the user through FUNC- 
TION statements. 

Initial conditions (other than zero) on integrators or transfer 
functions are specified with the INITIALIZE statement. 

There are currently five methods of numerical integration available 
in the MARSYAS system, Euler's l st -order, Butcher's St^-order, Sarafyan 
5th_ or der variable-step, 4 t ^ 1 -order Runga-Kutta and Adams-Bashforth 
predictor-corrector. These integration methods are discussed in Section 
IV B. The MARSYAS system automatically selects the Sarafyan variable- 
step method unless instructed otherwise by an INTEGRATE statement. If 
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the user chooses the Adams-Bashforth method, the MARSYAS system auto- 
matically selects a relative error of 0.002 unless otherwise specified. 

In addition, Runga-Kutta and Butcher's methods will be run using a step 
size of 0.01 unless otherwise specified. 

Conditions for temporarily or permanently halting the simulation 
are given by the STOP IF and TERMINATE IF statements, respectively. 

STOP IF is used to halt each simulation when using CHANGE statements. 

Constant gain elements (Constant Multipliers) can be changed into 
time-varying multipliers through the use of the VARY GAIN statement. 

If the user has described his system in block-diagram form in the 
Description Module, linear and/or nonlinear "loops" may be present. The 
absence of an integrator or transfer function in a closed path produces 
a loop. A loop is nonlinear if it contains at least one nonlinear ele- 
ment (such as a power function). Otherwise it is linear. In the MARSYAS 
system, loops are automatically solved using either a Newton-Raphson or 
a successive approximation technique and the elements within these loops 
are automatically listed for the user by the MARSYAS processor. The 
ESTIMATE statement can be used to assign initial values to the outputs 
of nonlinear elements within a nonlinear loop(s) should the user desire 
to do so, otherwise a value of zero is assumed. 

A summary of the statements which appear in the Simulation Module 
is shown in Figure 7 along with brief comments as to their function. 
Details of these statement operators are given in the following pages. 
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The order in which the statements appear in the M&RSYA.S deck is imma- 
terial except that SIMULATE must be the first statement of the Simula- 
tion Module and END the last, FUNCTION does not have to appear within 
the Simulation Module, but must appear before it is referenced in an 
EXCITE statement. It is recommended that all FUNCTION statements be 
placed at the very beginning of the MA.RSYA.S deck, after any DEVICE 
statements or just following the BEGIN statement. 
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OPERATOR 

CHANGE 

END 

ESTIMATE 

EXCITE 

FUNCTION 

INITIALIZE 

INTEGRATE 

PARAMETERS 

SIMULATE 

STOP IF 
TERMINATE IF 

VARY GAIN 

..L. 

A 

FIGURE 7. 


II C. 


FUNCTION 

Creates additional simulation runs with 
changes in the model. 

Marks the end of the Simulation Module. 

Specifies initial values for the outputs 
of nonlinear elements within a nonlinear 
loop. 

Specifies the excitation functions to 
be applied to the input terminals. 

Defines a non-standard excitation func- 
tion. 

Specifies initial conditions for inte- 
grators or transfer function elements. 

Specifies the integration method to be 
used. 

Specifies new parameters for a Descrip- 
tion Module element. 

Marks the appearance of the Simulation 
Module . 

Specifies the condition or conditions 
for stopping a simulation either tempo- 
rarily (STOP) or permanently (TERMINATE). 

\ 

i Defines time varying coefficients. 

Provides a comment for the user. 


STATEMENT OPERATORS USED IN THE SIMULATION MODULE 


79 



II C. SIMULATION MODULE 


CHANGE 


FORMAT - 


CHANGE: name of change; Statement operator : Statement format $ 

Note: For each change statement, a separate simulation is 

performed. 


COMMENTARY - 

The CHANGE statement is used when repetitive simulations are de- 
sired. For example, the user may wish to observe the behavior of his 
model when' certain parameters, excitation(s) , initial conditions or 
stopping conditions have been changed, A series of simulations can be 
set up in one run, thus saving valuable user and computer time. When 
using CHANGE, the user specifies an identifying name for the new simula- 
tion, the statement operator being changed and the statement format. 

The statement operators used with the CHANGE statement are PARAMETERS, 
EXCITE, INITIALIZE, STOP IF and TERMINATE IF. 

When changing Description Module ELEMENT parameters, the user must 

J 

specify the name of the MODEL in which the element occurs, the element 
name, the position number of the parameter in the parameter string and 
changed parameter values. 
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CHANGE (Cont'd.) 


EXAMPLE 

The following ELEMENTS statement appears in the Description 
Module model PSY 

ELEMENTS: TF, TF6 (2 , 1.3, 2.1, 7.4, 1.8, 3.4, 9.7)$ 

and the user desires to change the fifth value in the parameter string 
to 15.2. The CHANGE statement would be coded as follows: 

CHANGE: NEW TF6 ; PARAMETERS: PSY, TF6 (5 , 15.2) $ 

This causes the original ELEMENTS statement to be processed as 
if it had originally been coded as 

ELEMENTS: TF, TF6(2, 1.3, 2.1, 7.4, 15.2, 3.4, 9.7)$ 

If one of the parameters in an ELEMENTS statement is a MARSYAS 
name rather than a real number, the same technique as described above 
can be used to change its value. Consider the same ELEMENTS statement 
as shown in the above example except that the fifth element which we 
want to change has the name BILL as shown below: 

ELEMENTS: TF, TF6 (2 , 1.3, 2.1, 7.4, BILL, 1.8, 3.4, 9.7)$ 
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CHANGE (Cont'd.) 

Changing the value of BILL to 15 . 2 can be accomplished in the following 
manner : 


CHANGE: NEW TF6 ; PARAMETERS: PSY, TF6(5, 15.2)$ 

Since the name BILL is unique within the MODEL where it is used, 
the CHANGE statement can be shortened somewhat to take advantage of this 
fact by referencing the parameter name rather than the element name. The 
previous CHANGE statement could then be written as follows: 

CHANGE : NEW TF6 ; PARAMETERS : PSY , BILL ( 15 . 2 ) $ 

If several parameters of one element are to be altered using one 
CHANGE statement, the position numbers and new parameter values are 
listed in pairs, in any sequence, as shown below. 


EXAMPLE 

The following time varying coefficient element appears in the 
Description Module model TIM 

ELEMENTS: TV, TVl(3, 0.0, 0.0, 4.2, 1.6, 5.4, -3.2)$ 

and the user desires to change the last four parameters to the follow- 
ing values: 5.1, 1.93, 6.2, 1.1. The MARSYAS coding is as follows: 

CHANGE: RUN3 ; PARAMETERS: TIM, TV1(4, 5.1, 5, 1.93, 6, 5.4, 7, 1.1)$ 
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CHANGE (Cont'd.) 

An alternate method of coding the previous statement is as follows 

CHANGE: RUN3; PARAMETERS: TIM, TV1 (4,5.1) (5 ,1 .93) (6,5.4) (7,1 . 1) $ 

If an element contains only one parameter, the position number of 
the parameter is one. The coding for changing the value of a constant 
multiplier is shown in the example below. 


EXAMPLE 

If the constant multiplier element with name CM1 and parameter 
value 3.0 exists in the Description Module model ALPHA and the user 
desires to change the multiplier parameter value to 6.4, then the Simu- 
lation Module CHANGE statement with new simulation name NEW CM2 would 
read 


CHANGE: NEW CM2; PARAMETERS: ALPHA, CM2 (1, 6.4)$ 

If the user desires to change the Simulation Module EXCITE, 
INITIALIZE, STOP IF and TERMINATE IF statements, the form of the state- 
ment (s) is exactly as described elsewhere in this section of the Manual 
except that each statement is preceded by the word CHANGE and the title 
of the change. 
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CHANGE (Cont'd.) 


EXAMPLE 


CHANGE: name of change; EXCITE : input terminal name, new 

excitation function mnemonic$ 

CHANGE: name of change; INITIALIZE : model name, element name 

(new initial conditions) $ 

CHANGE: name of change; STOP IF : new logical express ion$ 

CHANGE: name of change; TERMINATE IF : new logical expression$ 

Except for STOP IF and TERMINATE IF, CHANGE statement altera- 
tions are permanent unless changed again in subsequent simulations. 

STOP IF statement changes are temporary and are valid only for the simu- 
lation specified by the given CHANGE statement. 

If there are "n" CHANGE statements in the Simulation Module, 
there will be "n+l" simulations executed by the MARSYAS processor. The 
additional simulation being that of the system as originally configured, 
as if no CHANGE statements were present. If the user desires to omit the 
original simulation and have the MARSYAS processor proceed to the simu- 
lation as specified by the first CHANGE statement, he must insert the 
following statement anywhere in the Simulation Module. 

CHANGE: (DELETE) $ 


84 



II C. SIMULATION MODULE 


CHANGE (Cont'd.) 

The simulations will then proceed, one by one, until all of the 
CHANGE statements have been executed o£ until a TERMINATE IF statement 
appears. 

When a TERMINATE IF statement is executed, the MARSYAS program is 
halted and no further CHANGE statements are executed. This is not the case 
with the STOP IF statement. For example, after changing one or more ele- 
ment parameters, initial conditions or excitation functions, the user may 
decide to examine the system outputs after a short period of simulation 
time, then proceed to the next CHANGE statement. The STOP IF statement is 
used for this purpose. 

EXAMPLE 


Let us assume that the user decides to change one of the system 
EXCITE functions to a ramp and also to change, the simulation time. The 
CHANGE statement could be coded as follows: 

CHANGE: RUN ONE ; EXCITE : INPUT6 , FRAMP (-2.0); STOP IF : TIME .GT. 10.0$ 

} 

I Note - Semicolon 

When there is more than one change specified in a particular 
CHANGE statement, each is separated by a semicolon as shown in the example 
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II C. SIMULATION MODULE 


END 


FORMAT - 


END$ 


COMMENTARY - 

The END statement is used to close the Simulation Module. The END 
statement must be the last one in the Simulation Module. 


EXAMPLE 


END$ 



II C. SIMULATION MODULE 


ESTIMATE 


FORMAT - 


ESTIMATE: model name, element name (output(s) estimate)$ 

COMMENTARY - 

The ESTIMATE statement is used to assign initial values to the 
outputs of non-linear elements which are part of the non-linear loop. 

When using the ESTIMATE statement, the user specifies the initial esti- 
mate of the output of a non-linear element which is part of a non-linear 
loop. The use of ESTIMATE is similar to that of INITIALIZE which speci- 
fies the initial conditions on integrators or transfer functions. The 
use of ESTIMATE is optional and is intended for use as an aid, if needed, 
in solving imbedded non-linear loops. 


EXAMPLE 

If a non-linear loop in model ALPHA contains a power function 
element whose name is PF1, and whose user specified estimate at time zero 
is 1.0, then the ESTIMATE statement would be 

ESTIMATE: ALPHA, PF1 (1.0)$ 
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II C. SIMULATION MODULE 


ESTIMATE (Cont'd.) 


Should a non-linear element contain two outputs, such as an out- 
put relay (with name R01) , then the ESTIMATE statement would be 


ESTIMATE: ALPHA, R01 (1.5, 3.2)$ 
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II C. SIMULATION MODULE 


EXCITE 

FORMAT - 

EXCITE: system input name, excitation function mnemonic (parameters) $ 

COMMENTARY - 

The EXCITE statement is used to specify an excitation function at 
each input terminal. The user must specify the excitation mnemonic and 
the parameters associated with the given excitation. A Table of Standard 
MARSYAS EXCITE statements appears at the back of this Manual. The Table 
lists the excitation function mnemonics and parameters, the type of ex- 
citation, the mathematical description of the excitation function, and 
a graph of the function versus time. The user is not restricted to the 
"standard" set. Non-standard excitations may be defined by the user 
through the use of the FUNCTION statement (see FUNCTION). When coding, 
non-standard excitations are treated in exactly the same manner as the 
standard MARSYAS set. 

EXAMPLE 

EXCITE: INPUT1 , FSIN (5.0, 2.0, 0.0)$ 

This statement causes the MARSYAS processor to excite terminal 
INPUT1 with the function 5 sin 2t. 
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II C c SIMULATION MODULE 

EXCITE (Cont'd.) 

The same excitation function can be used to excite several dif- 
ferent inputs as shown in the example below. 


EXAMPLE 


The system contains inputs INI, IN2 , IN3 which are to be excited 
with a step function of magnitude 3.2. The coding required to effect this 
is as follows: 

EXCITE: INI, IN2 , IN3, FSTEP(3.2)$ 

or 

EXCITE: INI, FSTEP(3.2)$ : IN2 , FSTEP(3.2)$ 

: IN3 , FSTEP(3.2)$ 
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II C. SIMULATION MODULE 


FUNCTION 


FORMAT - 

FUNCTION: a FORTRAN name, number of parameters $ 


COMMENTARY - 

The FUNCTION statement is used by the programmer to create a new 
excitation function, coded in FORTRAN, which is treated as though it 
were a MARSYAS Standard Excitation Function. This is done by specifying 
the new mnemonic as identified by the FORTRAN name, and the number of 
mathematical parameters required to compute the excitation function. A 
FORTRAN FUNCTION subprogram must immediately follow the MARSYAS FUNCTION 
statement as formatted above. The user-created FORTRAN FUNCTION subpro- 
gram then computes the output of the new excitation function from the 
parameter values specified in the EXCITE statement. (See EXCITE.) 

The subprogram name referenced in the FORTRAN FUNCTION statement 
is the same FORTRAN name referenced in the MARSYAS FUNCTION statement. 

This is the excitation function mnemonic assigned in the Simulation Module 
EXCITE statement. The arguments in parenthesis which follow the FORTRAN 
FUNCTION name are arbitrary real FORTRAN names for the parameters re- 
quired to compute the excitation function in the subprogram. The last argu- 
ment in the parameter list must be a symbol for time. The value of time 
must not be altered or otherwise manipulated within the 1 subprogram. The 
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IT C. SIMULATION MODULE 

FUNCTION (Cont'd.) 

excitation name referenced in the MARSYAS EXCITE statement is a MARSYAS 
name and, hence, is not restricted to six characters as are FORTRAN 
names . 

It is permissible to place the FUNCTION statement (s) and associated 
subroutine (s) anywhere within the MARSYAS program after the BEGIN state- 
ment and before the Post Processing Module provided they appear prior to 
being referenced in EXCITE statements. It is recommended, however, that 
all FUNCTION statements be placed at the beginning of the deck, after or 
shortly following, the BEGIN statement to isolate the source FORTRAN 
code from the MARSYAS language statements. 


EXAMPLE 

Shown below is the FORTRAN coding for the truncated ramp excitation 
function shown. Following it is the MARSYAS coding referencing this 
FUNCTION. 
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Coding 
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FUNCTION: RAMP, 3 $ 

FUNCTION RAMP (A, B, T) 

IF (( T .GE. 0.0) o AND. (T .LT. B)) 
IF (T .GE. B) X = 0.0 
RAMP = X 
RETURN 


SIMULATE: EP0CH$ 

• • • • • 

EXCITE: INPUT6 , RAMP (2.0, 2.0)$ 

• • • • • 

END$ 
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X = A/B * T 



II C. SIMULATION MODULE 


* 


• INITIALIZE 

FORMAT - 

INITIALIZE: model name, element name (initial conditions) $ 

COMMENTARY - 

The INITIALIZE statement is used to assign initial conditions to 
the integrators and transfer functions defined in the MARSYAS program. 
When using the block-diagram approach, the user must specify the name 
of the model in which the integrator or transfer function appears along 
with its element name and initial conditions. If specifying the initial 

■hVi 

conditions of an n--order transfer function, n-initial conditions must 
be specified in the ascending order of the derivatives as shown in 
EXAMPLE (1). If the INITIALIZE statement is omitted, all initial condi- 
tions are defaulted to zero, automatically. 

When using the EQUATION mode of solution as defined in the Descrip- 
tion Module, the initial conditions are specified by stating the model 
name followed by each derivative and its initial condition. If one of 
the derivatives is omitted, then its initial condition is automatically 
defaulted to zero. See EXAMPLE (2). 

When using a model which has a MATH REFERENCE transfer function, 
the initial condition on the transfer function is specified as shown in 
EXAMPLE (3). 
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II C. SIMULATION MODULE 


INITIALIZE (Cont'd.) 

EXAMPLE (1) 

If model ALPHA contains an integrator with name IN].,, initial condition 
of 3.6, and second order transfer function with name TF6, where TF6 is 

= S-Li with y (0) = 1.2 and y' (0) = 4.4 

s 2 + 2s + 1 

then the initial conditions on each are specified as follows: 

INITIALIZE: ALPHA, INI (3.6)$ 

: ALPHA, TF6 (1.2, 4.4)$ 

EXAMPLE (2) 

Suppose in model BETA, which uses EQUATION, the variable X as well 
as the derivatives X"", X"' , X", and X* appear. Then the non-zero ini- 
tial conditions would be expressed as follows: 

INITIALIZE: BETA, X m (1.6), X" (2.2)$ 

The initial condition on X' would automatically be set to zero 
since it has been omitted from the INITIALIZE statement. 
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II C. SIMULATION MODULE 


INITIALIZE (Cont* d. ) 

EXAMPLE (3) 

Suppose we wish to code the following expression in EQUATION form 
in a Description Module model whose name is PHY. 

Y(s) = 1 

X(s) s + 1 

Using the MARSYAS Transfer Function element as a MATH REFERENCE, 
the above could be written in EQUATION form as follows: 

MODEL: PHY, EQUATION$ 

• • • • 

Y = TF (1, 0, 1, 1, 1, X) $ 

• • • • 

END$ 

Then, in the Simulation Module, the INITIALIZE statement would read 
INITIALIZE: PHY, Y (1.6)$ 
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II C. SIMULATION MODULE 

INTEGRATE 

FORMAT - 

INTEGRATE: mnemonic, code word, real number $ 

COMMENTARY - 

The INTEGRATE statement (optional) is used to specify the integra- 
tion algorithm to be used in the problem solution. Five methods of inte- 
gration are presently available in the MARSYAS system, Euler's (l st -order) 
Butcher's (5 th -order), Sarafyan (5 th -order) variable-step, Runga-Kutta 
(4th_ or der) and Adams-Bashforth predictor-corrector. Other numerical 
integration methods soon will be added to the MARSYAS language allowing 
the user greater flexibility in solving his problem. Only one INTEGRATE 
statement can appear in a MARSYAS deck. When choosing the method of inte- 
gration desired, the user specifies its mnemonic, the code word TIMESTEP, 
and integration step size. The mnemonics for each of the integration 
methods currently available are as follows: Euler's (EU) , Butcher's (BU) , 

Sarafyan (SA) , Runga-Kutta (RK) , Adams-Bashforth (AB) . The timestep may 
be specified only for Euler, Butcher and Runga-Kutta methods. If the 
timestep is omitted, the simulations will be performed with an integra- 
tion step size of 0.01. 

If the user omits the INTEGRATE statement, the simulation will be 
run using the. Sarafyan method. If the user chooses the Adams-Bashforth 
method, the MARSYAS system automatically selects a relative error of 0.002 


97 



II C. SIMULATION MODULE 


INTEGRATE (Continued) 


An absolute error or other relative error may be selected if the code 
word ABSERR or RELERR, respectively, is used in an INTEGRATE statement, 
followed by the user specified error. 

If the simulation problem contains a Sample and Hold element(s), 
Butcher's method or the Runga-Kutta method of integration should be used, 
with the TIMESTEP smaller than the smallest Sample and Hold interval. 


EXAMPLES 


a) 

b) 

c) 

d) 

e) 


INTEGRATE: 

RK, 

TIMESTEP, 0.01$ 

INTEGRATE: 

AB, 

RELERR, 0.05$ 

INTEGRATE: 

AB, 

ABSERR, 0.001$ 

INTEGRATE: 

BU, 

TIMESTEP , 0.01$ 

INTEGRATE: 

EU, 

TIMESTEP, 0.01$ 


Euler's method is implemented automatically whenever the MARSYAS 
processor detects the presence of a discontinuity. Thus, if a disconti- 
nuity occurs at time, t^, the integration method will change to the 
Euler method at time t^-€ and back to the method of integration used before 
the discontinuity was encountered, at time t^+€. The reason for the 
Euler-interrupt is that the lower order method gives better results at 
discontinuities than the higher order methods. See Section IV B. 

To disable the Euler-interrupt at discontinuities so that the chosen 
integration scheme continues without interruption, simply insert the 
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II C. SIMULATION MODULE 


INTEGRAT E (Continued) 


word XEULER after the timestep or error value in the INTEGRATE statement 
as follows: 


INTEGRATE: RK, TIMESTEP, 0.001, XEULER$ 
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PARAMETERS 


FORMAT - 

PARAMETERS: model name, element name (parameter number, numerical 

value) $ 
or 

PARAMETERS: model name, parameter name (numerical value) $ 

COMMENTARY - 

The PARAMETERS statement is used to assign numerical values to 
Description Module element parameters in much the same way as the 
Description Module PARAMETERIZE statement (Section II A.). The only 
difference between the two operators is that the model name must be 
specified in the PARAMETERS statement, whereas it is not required in the 
PARAMETERIZE statement. 

The user might choose to use this operator when altering parameters 
in several different Description Modules thus eliminating the need for 
changing his original coding. 

EXAMPLE 

PARAMETERS: PSY, CM1 (1, 2.2)$ 

See PARAMETERIZE and CHANGE for additional comments and examples. 
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II C. SIMULATION MODULE 


SIMULATE 


FORMAT - 


SIMULATE: System model name$ 


COMMENTARY - 


The SIMULATE statement is used to identify the system being 
simulated. The system model name is the main model name specified in 
the Description Module. Since there is only one system model, any 
other models appearing in the Description Module describe the structure of 
submodels appearing in the main system model. 

SIMULATE must be the first and END the last statement of this 

module. 

EXAMPLE 

SIMULATE: SPACE SHUTTLE ENGINE$ 
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II C. SIMULATION MODULE 


STOP IF 


FORMAT - 


STOP IF: logical expression $ 

COMMENTARY - 

The STOP IF statement specifies the condition(s) for temporarily 
halting a simulation. When more than one simulation is to be executed 
(via CHANGE statements), a halt caused by STOP IF allows the execution 
of the next simulation to proceed. This is contrasted with a halt 
caused by TERMINATE IF which permanently stops the simulation and pre- 
vents the execution of any additional simulations. There is no differ- 
ence between STOP IF and TERMINATE IF if no CHANGE statements are present. 
The construction of the logical expression part of the STOP IF statement 
is described under TERMINATE IF. There may be as many STOP IF statements 
as the user may require to control the simulation. 

\ 

i 

EXAMPLE 

STOP IF: TIME .GT. 20.0$ 

STOP IF: (0UTPUT6 .LE. 0UTPUT3)$ 

STOP IF: (INPUT1 .AND. INPUT2) .GE. 60.0$ 
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II C. SIMULATION MODULE 


TERMINATE IF 


FORMAT - 


TERMINATE IF: logical expression $ 

COMMENTARY - 

The TERMINATE IF statement specifies the condition (s) for perma- 
nently halting a simulation. The halting conditions are expressed 
using standard FORTRAN logical expressions. The logical operators AND, 
OR, NOT and relational operators EQ, NE, GT, GE, LE and LT are used to 
form logical expressions involving TIME, and/or model INPUT or OUTPUT 
names. A halt caused by TERMINATE IF halts the simulation being exe- 
cuted and no additional simulations will be processed even though 
CHANGE statements may be present. (See STOP IF.) There may be as many 

TERMINATE IF statements as the user may require to control the Simula- 

/ 

tions. 

\ 

\ 

EXAMPLE 

TERMINATE IF: TIME .GT. 5.0$ 

TERMINATE IF: (OUTPUT1 .OR. 0UTPUT2 ) .GT. 3.2$ 

TERMINATE IF: (INPUT1 .AND. INPUT2) . EQ. 60.0$ 

Testing on equality as shown above should be avoided because the 
equality may never be exactly satisfied. ; 
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VARY GAIN 


FORMAT - 


VARY GAIN: 


model name, 


constant multiplier 
element name 


( algebraic expression 
or numerical list 


COMMENTARY - 

The VARY GAIN statement is used to change Constant Multiplier elements 
into Time Varying multipliers. This feature permits the user to solve 
different equations with time varying coefficients or solve models which 
contain time varying multipliers. Recall that Constant Multiplier (CM) 
elements are assigned names in a DESCRIPTION Module using an ELEMENTS state- 
ment. These CM element names along with the model name of the DESCRIPTION 
Module in which they appear are used in the VARY GAIN statement above. 

EXAMPLE 

Suppose in model BETA, with Constant Multiplier element names 
MULTI (36. 3), MULT2(0.05) and MULT3(-3.4), the user desires to change 
MULT2 to 3t^ +5.2 cos t. Then the VARY GAIN statement would be expressed 
as follows: 

VARY GAIN: BETA, MULT2('3.0 * TIME ** 2 + 5.2 * COS (TIME) ')$ 

Note the quote marks ' ' which must enclose the entire expression. 
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II C. SIMULATION MODULE 


VARYjGAIN (Continued) 

Should the time-varying nature of the function be expressed as a 
linear piecewise continuous curve, then the coordinates of the curve at 
the discontinuities are sufficient to define the function as shown below 



which is represented by the numerical list 
( n, t 1 ,y 1 , t 2 ,y 2 , t^) 

where the first number in the list, n, is an integer specifying the num- 
ber points defining the function (^100 points), and the numbers following 
n are the coordinates of the points. Note that the numerical list is 
not enclosed by quote marks. 

EXAMPLE 

In the previous example, let MULT2 be a piecewise continuous func- 
tion similar to the plot above. Then the VARY GAIN statement would read 
VARY GAIN: BETA, MULT 2 (4, 0.0, 1.0, 1.2, 2.2, 2.8, 1.0, 4.3, 2.4)$ 
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D. Post Processing Module 

The function of the Post Processing Module, which always follows 
the Simulation Module, is to specify the output format for the MARSYAS 
processor. In addition, the Post Processing Module contains a Fast 
Fourier Transform (FFT) processor which is explained later in this sec- 
tion. A MARSYAS program can contain only one Post Processing Module. 

The following is a brief description of the terminology which the 
user will encounter in reading this section of the Manual. 

The PRINT statement is used to specify which output variables are 
to be printed while the SAMPLE statement indicates the print interval. 

If all or part of the output is to be plotted, a PLOT statement is re- 
quired. Should the frequency response at any output terminal be desired, 
the FOURIER statement is used. 

A summary list of the statement operators used in this section 
appears in Figure 8 along with brief comments as to their function. The 
order in which the statements appear within the Post Processing Module 
is immaterial except that END must be the last statement. 

A more detailed description of these operators as well as examples 
as to their use appears in alphabetical order on the following pages. 
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OPERATOR 


FUNCTION 


END 

Marks the end of the Post Processing 
Module. 

FOURIER 

Calculates the frequency response 
of the system inputs or outputs 
specified. 

PLOT 

Defines the variables to be plotted 
and the way they are to be plotted. 

PRINT 

Specifies the variables to be printed. 

SAMPLE 

Defines the sampling rate at which the 
variables are printed. 

* 

Provides comment for the user. 


FIGURE 8. STATEMENT OPERATORS USED IN THE POST PROCESSING MODULE 
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END 


FORMAT - 


END$ 


COMMENTARY - 


This should be the last statement in the post-processing module 



II D. POST-PROCESSING MODULE 


* 


FOURIER 


FORMAT - 

FOURIER: name 1 (period length), name 2 (period length), etc. $ 

COMMENTARY - 

The FOURIER statement provides the user with a frequency response 
of the system inputs or outputs specified. The frequency response con- 
sists of a Fourier Transform magnitude and phase spectrum and power 
spectral density (PSD). The waveform period length must be specified 
by the user and is indicated by a real number in parenthesis after the 
input or output terminal name. A maximum of 2048 samples are available 
for use in obtaining the frequency response of each terminal. Since 
Runga-Kutta is the only integration scheme permitted when using FOURIER, 
the user must coordinate his integration step size with the period length 
to take advantage of the maximum number of sample points which the MARSYAS 
processor is capable of handling. In addition to tabulation data output, 
the user is provided with linear plots of PSD, magnitude and phase angle 
versus frequency. 

EXAMPLE 

The user has decided to obtain the frequency spectrum of his sys- 
tem at output LAMBDA. After having examined the time response at LAMBDA, 
he decides that a period length of 16.0 seconds is required to define 
his waveform. In order to utilize all 2048 permissible sample points. 
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• FOURIER (Cont'd.) 

the user computes his Runga-Kutta step size to be 

16.0 s econds -*■ 2048 ^ 7-7 = 0.0078125 seconds/ sample point, 

period period r r 

This is the integration TIMESTEP which must be used as part of the 
INTEGRATE statement in the Simulation Module. The FOURIER statement is 
coded as follows: 

FOURIER: LAMBDA (16.0)$ 

It is not necessary for the user to use all 2048 sample points in 
order to obtain the frequency response. However, the use of the maxi- 
mum number of points (2048) will always give more accurate results than 
when using fewer points. 
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PLOT 


FORMAT - 

PLOT: gn d -type (TIMEl,TIME2,CODE) , N(XL,Y1, X2,Y2, X3,Y3)$ 

COMMENTARY - 

The PLOT statement is used to obtain automatically scaled plots 
of the INPUTS and/or OUTPUTS of the system being simulated. The grid 
type is specified by the mnemonics LINEAR, SEMILG, LOGLOG and LOGSEM. 
Figure 9 summarizes the functions of the grid type mnemonics. TIME1 
and TIME2 are the plot start and stop reference times. CODE is an inte 
ger, 1 or 2, indicating the nature of the independent variable (s) being 
plotted. If CODE is 1, the independent variable being plotted is TIME; 
and, if CODE is 2, it is not TIME. Code equal to one may be implied, 
if desired, by closing the parenthesis after TIME2. 

The integer N is a positive or negative number which designates 
the number of frames of graph paper to be used when plotting the data. 

N can be greater than one only when TIME is the independent variable 
being plotted. If N is one, both N and the parenthesis may be omitted. 

A frame is a unit of output paper approximately 7.5 inches square. 
Two frames would be 7.5 inches x 15 inches; three frames would be 
7.5 inches x 22.5 inches; etc. If the user decides to spread out the 
TIME axis, he can do so by specifying as many frames, N, as needed to 
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PLOT (Cant'd.) 


GRID TYPE 
MNEMONIC 

ABSCISSA 

ORDINATE 

LINEAR 

1 ine ar 

linear 

SEMILG 

linear 

logarithmic 

LOGLOG 

logarithmic 

logarithmic 

LOGSEM 

logarithmic 

linear 


FIGURE 9. AXIS DESIGNATION FOR VARIOUS PLOT MNEMONICS 
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PLOT (Cont’d. ) 

Normally, there appears one set of arid lines per frame of graph 
paper. There may be times, however, when the user decides to place more 
than one grid, one above the other, on a frame of graph paper. A maxi- 
mum of three grids can be placed on one frame of graph paper. A nega- 
tive sign preceding the integer N denotes that multiple grids are de- 
sired. The number of coordinate names in parenthesis following N denotes 
the number of grids needed. Since a maximum of three grids are permitted 
on a frame of graph paper, up to three sets of coordinate names are per- 
mitted inside the parenthesis. A maximum of three plots can be placed 
on one grid if multiple plots are desired. A positive sign or lack of 
sign for the integer N denotes that multiple plots are required with the 
number of sets of coordinate names in parentheses indicating the number 
of plots needed. 

XI and YI are the user-given MARSYAS names of the system INPUTS 
or OUTPUTS which are to be plotted. XI is the independent variable (the 
abscissa), TIME or other system OUTPUT or INPUT name, and YI the depend- 
ent variable (the ordinate). A maximum of three sets of variable names 
may appear within the parenthesis of a given PLOT statement if paren- 
theses are used. See EXAMPLE (1). As many PLOT statements as required 
may be used within the Simulation Module so long as the total number of 
graphs processed does not exceed one hundred. 
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PLOT (Cont'd.) 

If a simple plot of one variable against another on one frame of 
graph paper is desired, the integer N and the parenthesis may be omitted 
from the plot statementc Then, up to one hundred coordinate names may 
be plotted using a single PLOT statement, as 

PLOT: grid type (TIMEl ,TIME2 ,CODE) , Xl,Yl, X2,Y2, . .., X100,Y100$ 

mnemonic 

(See EXAMPLE (2).) 

Additionally, if the independent variable being plotted is time 
(CODE K 1), the abscissa name, TIME, need not be specif ied, thus further 
simplifying the PLOT statement. The PLOT statement would then become 

PLOT: grid type (TIMEl, TIME2 ) , Yl, Y2, Y3, ..., Y100$ 

mnemonic 

(See EXAMPLE (3).) 

EXAMPLE (1) 

The user decides to cross plot three OUTPUT variables with names 
VI, V2, V3 on a single sheet (one frame) of linear graph paper using three 
separate grids. The reference time interval will be from 0 to 36 seconds. 
The variables are to be plotted in the following manner: V2 vs. VI, 

V3 vs. V2 and VI vs. V3. 

The PLOT statement is coded as follows: 

PLOT : LINEAR (0.0,36.0,2), -1 (VI, V2, V2,V3, V3,V1)$ 
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II D. POST- PROCESSING MODULE 


PLOT (Cont'd.) 

Note that N s -1 indicates that one frame of graph paper is to be 
used. The negative sign indicates multiple grids are needed. The three 
coordinates inside the parenthesis indicate that three grids are required. 
The execution of the above statement in an actual MARSYAS program is 
shown in Figure 10 . 

EXAMPLE (2) 

The user decides to cross plot OUTPUTS V2 and VI with Vl the 
abscissa on a single frame of linear graph paper in the time interval 
0 to 36 seconds. Since N=l, both N and the parenthesis may be omitted 
from the PLOT statement as shown below. 

PLOT : LINEAR (0.0,36.0,2), VI, V2$ 

The execution of the above statement in an actual MARSYAS program 
is shown in Figure 11. 

EXAMPLE (3) 

A plot of system OUTPUT V3 against TIME is desired, on one frame 
of graph paper with a linear grid. The time interval is 0 to 36 seconds. 
Since N=1 and TIME is the independent variable, the abbreviated form of 
the PLOT statement is used. Figure 12 is an actual MARSYAS execution of 
the following statement: 

PLOT : LINEAR(0.0, 36.0, 1), V3 $ 
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PLOT (CON’T) 



FIGURE 12 PLOT OF VARIABLE .AGAINST TIME 
PLOT: LI^EAR(0„0,36c0,T),V3 $ 
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PRINT 


FORMAT - 

(1) PRINT: name 1, name 2, etc. $ 

or 

(2) PRINT: (output header title) $ 

COMMENTARY - 

The PRINT statement specifies the names of the system input and/or 
output data to be printed. FORMAT (1) above yields the standard MARSYAS 
output print listing of six columns of data per page. The first printed 
column is always the independent variable TIME. FORMAT (2) allows the 
user to specify a title heading which will appear at the top of the out- 
put printed page. 

The appropriate SAMPLE statement must be used in conjunction with 
PRINT to indicate the print step desired (see SAMPLE). There is a limit 
of one hundred input and output names that may be included in a given 
PRINT statement. The MARSYAS processor will automatically print the 
indicated variables in groups of six until the output list is completed. 

EXAMPLE 

(1) Assume the usep's system has three inputs, INPUT1, INPUT2, 
and INPUT3, and two outputs, 0UTPUT1 and 0UTPUT2, which he wishes to be 
printed out, along with the header title FIRST TEST CASE. He decides to 
examine his input-output data at every fifth integration step. The 
PRINT and associated SAMPLE statements would be coded as follows: 
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PRINT (Cont'd.) 


PRINT: (FIRST TEST CASE)$ 

PRINT: INPUT1 , INPUT2, INPUT3, 

OUTPUT1 , OUTPUl'2 $ 
SAMPLE: STEP, 5$ 


When using Runga-Kutta, if instead of examining the output at every 
fifth integration step the user decides to examine the data at every 
0.10 seconds, the PRINT statement would remain unchanged and the SAMPLE 
statement becomes 


SAMPLE: TIME, 0.10$ 


I 
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II D„ POST-PROCESSING MODULE 

SAMPLE 


FORMAT - 

(1) SAMPLE: STEP, integer number $ 

or 

(2) SAMPLE: TIME, real number $ 

COMMENTARY - 

The SAMPLE statement is used in conjunction with the PRINT state- 
ment to specify the output printing interval. The print interval may 
be expressed in terms of the independent variable TIME or the integra- 
tion STEP size. If TIME is selected, the step size is taken in approxi 
mately equal increments of simulation time as specified by the real 
number in Format (2) above. 


EXAMPLE 


SAMPLE: STEP, 5$ 

or 

\ 

SAMPLE: TIME, .05$ 

In the first example, the MARSYAS processor will list the output(s) 
at every fifth integration step, whereas the second will list the out- 
put (s) in time increments of 0.05 units. 
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SECTION III 

SOFTWARE STRUCTURE OF MARSYAS 
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III. SOFTWARE STRUCTURE OF MARSYAS 

A. Introduction 

The prime objective of the MARSYAS software is to transform a MARSYAS 
program, which describes a model and specifies the simulation, into a 
FORTRAN program that contains the arrays and subroutines for the numeri- 
cal solution of the various matrix equations. The MARSYAS software is 
thus a precompiler which compiles the MARSYAS language statements into a 
set of FORTRAN programs called the Object Program which is then executed 
by the computer. 

B. Overview of MARSYAS Software System 

MARSYAS was originally designed for use on a time-sharing machine 
such as the Univac 1108 computer. The software permits several users to 
access MARSYAS simultaneously from remote stations. An overview of the 
MARSYAS software system is shown in Figure 13. The user's block diagram 
and/or equations are coded in a MARSYAS program consisting of Description, 
Modification, Simulation and Post Processing Module statements which are 
punched on cards and fed to the computer. All Program Modules with the 
exception of the Post Processing Module feed into the Simulation Program 
Module which ultimately generates the FORTRAN Object Program. This, in 
conjunction with the Library of Standard Elements and Excitation Functions, 
is compiled with the FORTRAN Object Program and executed. 
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ill B. OVERVIEW OF MARSYAS SOFTWARE SYSTEM 


MARSYAS- PROGRAM 




EXECUTION of FORTRAN 
OBJECT PROGRAM 


SIMULATION 3'JN i 

I 


PRINT -O'JT, PLOTTING 


FIGURE 13. OVERVIEW OF MARSYAS SOFTWARE SYSTEM 


















Ill C. LIMITATIONS OF MARSYAS 


C. Limitations of MARSYAS 

Within the complete MARSYAS program, there is one main MODEL and 
up to 99 SUBMODELS permitted. Each SUBMODEL, of course, has all the 
characteristics of a MODEL. Models and Submodels are each limited to 
a maximum of 50 inputs, 50 outputs, 300 elements with restrictions as 
noted below and 500 connections. 

For the main model and all its submodels, the following restric- 
tions apply: 

There may be a maximum of : 

350 constant multipliers 

150 transfer functions, including integrators 

If is the order of each transfer function, M 
is the total number of transfer functions and 
Lf = Ni + 1, then 

M 

12 L i ^ 400 

i=l 

v 

400 nonlinear elements 
4000 parameters 

600 inputs to nonlinear elements 

400 outputs from nonlinear elements 
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Ill c. ’limitations of marsyas 

For the Simulation Module, the following restrictions apply 
There may be a maximum of: 

125 EXCITE statements 

600 STOP IF and TERMINATE IF statement words 
30 simulation CHANGE cycles 
400 parameter changes per cycle 

A "word” is the number of characters divided by six. 

The following restrictions apply to the Simulation Module: 
There may be a maximum of: 

100 plot terminals 

100 graphs 

100 print terminals 


350 Hollerith words for PRINT title statement 
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MATHEMATICAL STRUCTURE OF MARSYAS 



IV. MATHEMATICAL STRUCTURE OF MARSYAS 

A. The State Space Approach 

MARSYAS has a basic mathematical structure which incorporates the 
latest and most up-to-date mathematics associated with modern control 
theory. The state space approach is used throughout, and directed graph 
theory is used to detect and unwind algebraic loops automatically. These 
methods have considerable advantages over older, classical methods. 

There are no known methods that will solve every set of nonlinear 
algebraic equations. Therefore, while many nonlinear loops can be 
solved by MARSYAS, some will arise which will not converge. However, 
most all cases of linear loops can be handled. An overview of the mathe- 
matical process which is the foundation of MARSYAS appears in Figure 14. 

B. Numerical Integration Techniques 

The FORTRAN Object Program generated by the MARSYAS processor uses 
precoded subroutines to control the numerical integration and detect 
discontinuities which may exist in the model. Since serious errors may 
arise when integrating through discontinuities, the MARSYAS system auto- 
matically switches integration methods when a discontinuity is detected. 
It integrates through the discontinuity using Euler's method. Euler's 
method, being of first order, will yield more accurate results at dis- 
continuities than higher order methods. Once past the discontinuity, 
the original integration scheme is resumed. 
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IVA. THE STATE SPACE APPROACH 
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FIGURE 14, OVERVIEW OF MATHEMATICAL PROCESS 















IV B. 'NUMERICAL INTEGRATION TECHNIQUES 


1 . Runga-Kutta Method 

The Runga-Kutta-Gill numerical integration scheme used by 
MARSYAS is the fourth-order single step method which approximates the 
solution of the state equation X(t) = F(X(t), t) at t = (n+l)h by 

X n+1 = X n + I (K 1 + 2 (1 ' ) ^ + 2 (1 + s/TS) ^ 

where 

K 1 = »C*n. tn> 

K 2 = F(X n + K x h/2, t n + h/2) 

K 3 = FCXn - (% - V3 ) Kjh + (1 - V75 ) K,h, t n + h/2) 

K4 = FCXj, - y/~J K 2 h + (1 + VTf) K 3 h, t n + h) 

and h is the step size specified by the user. 

2. Adams-Bashforth Method 

MARSYAS actually uses the Adams-Bashforth predictor with the 
Adams-Moulton corrector in solving the state equations. These multi-s 
equations approximate the solution of the state equation as follows: 

Adams-Bashforth predictor: 

*n+l = *n + h (55 *n " 59 *n-l + 37 *n-2 ‘ 9 *n-3> 



IV Bo NUMERICAL INTEGRATION TECHNIQUES 


Adams-Moulton corrector: 

x n+l = x n + (^n+l + 19x n “ 5x n-l + x n-2> 

These predictor-corrector equations have a definite advantage 
with respect to solution time over the Runga-Kutta method in that at 
most two derivatives must be calculated at each time step whereas the 
Runga-Kutta method requires four derivatives. The automatic error 
estimation procedure inherent in the predictor-corrector scheme allows 
the predictor to select the optimum step length that satisfies the built- 
in error criterion, thus further speeding up the solution time. 

3 . Euler's Method 

Euler's method is implemented automatically whenever the MARSYAS 
processor detects the presence of a discontinuity. Euler's is essentially 
a first-order single step method of the form 

X(n+1) = X(n) + X(n) h 

The integration step size h is measured from the end of the 
last full integration step before the discontinuity and extends past the 
discontinuity as shown in the sketch below. 
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IV B. NUMERICAL INTEGRATION TECHNIQUES 


Thus, if a discontinuity occurs at time -t=b, then the integra- 
tion method will change to the Euler method at t e t and back again at 

t=t , . to the method of integration used before the discontinuity was 
n+1 

encountered. 

If the user has several discontinuities in his system, it is 
recommended that Butcher's method or the Runga-Kutta method of integra- 
tion be specified (using the INTEGRATE statement) since the solution time 
will often be less than when using Adams-Bashf orth. 

Should the need arise, the user can disable the Euler-interrupt 
so that the selected integration scheme can proceed uninterrupted. Addi- 
tionally Euler's method can be used for the entire simulation if desired. 
See INTEGRATE statement in Section II C. 

4. Butcher's Method 

The Butcher's numerical integration method use in MARSYAS is a 
fifth-order, six-stage Runga-Kutta scheme which approximates the solution 
of the state equations as follows: 

x .. ■ x + -r—r (7K + 32K_ + 12K. + 32K C + 7K,) 

n+1 n 90.0 1 2 4 5 6 

where 

I 

K, « F (x ,t ) 

1 n’ n 

K„ <= F(x + KMi/4, t + h/4) 

Zonin' 

K„ = F (x + k h/8 + K.h/8, t + h/4) 

3 n 1 2 

K 4 = F (x n - K 2 h/2 + K 3 h, t + h/2) 
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IV B. NUMERICAL INTEGRATION TECHNIQUES 




K 5 = F(x n + 3/16K 1 h + 9/16K 4 h, t + 3/4h) 

K, = F(x - 3/7K,h + 2/7K 0 h + 12/7K„h - 12/7K h + 8/7K h, fc + h) 

6 n 1 2 3 h -> 

and h is the step size specified by the user, 

5. Sarafyan Method 

The Sarafyan variable step method is essentially an embedding 
of Runga-Kutta formulas to achieve step size control. The Sarafyan 
method used in MARSYAS is essentially a fifth-order, six-stage Runga-Kutta 
formula (Butcher's) with an imbedded second-order Runga-Kutta formula used 
to determine the appropriate step size. The fifth-order, six-stage formula 
is the same as that used in the Butcher method above and the second-order 
formula is given by 


x 


n+1 


x + h(-K + 2K ) 
n J. z 


where 

K n = F (x , t ) 

Inn 

K„ = F (x + K.h/4, t + h/4) 

2 n I n 

and h is the step size. 

The higher and lower order Runga-Kutta schemes are compared 
so that an estimate of the accuracy can be made. Then a judgment is 
made on the current value of h being used. If h is too large, it is 
halved. If too small, it is doubled. If just right, it is left as it is. 

Computation accuracy using the Sarafyan method is high and 
computation speed fast. Depending upon problem being simulated, computa- 
tion speed can be increased significantly over fixed-step methods. 
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Other numerical integration methods will soon be added to the 
MARSYAS language allowing the user greater flexibility in solving his 
problems , 
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IV C. SOLUTION SCHEMES FOR DIFFERENTIAL EQUATIONS 

C. Solution Schemes for Differential Equations 

When using the EQUATION option in solving a set of ordinary dif- 
ferential equations, the user must check to see that a solution exists. 

In other words, it must be possible to construct a "solution scheme" for 
the system of equations. A solution scheme is an assignment of variables 
to the equations such that the following three conditions hold: 

1. There must be one independent equation for each unknown variable. 

2. The highest order of each variable derivative in the set of 
equations can be "assigned" to one particular equation of the 
set. No two variable derivatives can be assigned to the same 
equation. 

3. The order of the variable derivative assigned to a particular 
equation must be the highest order for that variable as it 
appears in the system of equations. 

EXAMPLE (1) 

Determine if a solution scheme exists for the following set of 

I 

differential equations: 

1) aX + bYX + Z = 0 

2) Y + Z = 2XY 

3) Z + X = 0 

In the above set of equations in three unknowns, X, Y and Z, the 
three highest order derivatives are X, Y, and *Z. They can be assigned 
to equations 1), 2) and 3), respectively. Hence a solution scheme exists 
for the system of equations. 
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> 

EXAMPLE (2) 

Determine if a solution scheme exists for the following set of 
differential equations: 

1) X + Y + XY = SIN (X - Y) 

2) XY - XY = 0 

These are two equations in two unknowns, X and Y, whose highest 

• • • •• • • 

order derivatives are X and Y. The derivative X can be assigned to 
equation 1) and Y can be assigned to equation 2). Hence, a solution 
scheme exists for this set of equations. 

EXAMPLE (3) 

Determine if a solution scheme exists for the following set of 
differential equations: 

1) X + Y + XY = TAN (3X 2 + 2Y) 

2) XY - X = 0 

The above two equations in two unknowns X and Y have highest 

• • • 

order derivatives X and Y. Since we cannot assign two variable deriva- 
tives to the same equation, no solution scheme exists for this set of 
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V. OPERATION OF THE MARSYAS SYSTEM 

A. Deck Setup for MARSYAS Operation on the Univac 1108 Computer Under 

EXEC VIII 

The following computer control cards are required to run a MARSYAS 
program on the Univac 1108. 

@RUN, 

@ASG , X MARSYAS'-MARSYAS 

@ASG ,T 1 ,F 

0DATA, I 1 

MARSYAS DECK 

@END 

@ADD MARSYAS “MARSYAS . 

@FIN 

B. MARSYAS Diagnostics 

There are two distinct types of error diagnostics in the MARSYAS 
system. Statement Error Messages and Module Error Messages. The first 
type deals with errors encountered in coding a given statement. When 
an error of this type is detected, the processing of the MARSYAS state- 
ment is terminated with an error message. That portion of the statement 
in error is indicated with a series of backward slashes ( \ ) under the 
printed out statement. An error diagnostic is then printed out to aid 
the user in correcting the problem. Since the MARSYAS processing termi- 
nates after processing a group of errors, the user should recheck his 
program to be sure no more errors are present in the non-processed por- 
tion of his program. 
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The second type of error diagnostic will appear when the user vio- 
lates the structural pattern of MARSYAS. Examples of this are if the 
user should inadvertently place the Post Processing Module before the 
Simulation Module, or if he forgets to place an END$ statement at the 
end of a Module and so on. If a system is incompletely defined by 
leaving out elements, connections, inputs, outputs, excitations, and so 
on, the errors will be grouped at the end of the appropriate Module 
under the heading "MODULE ERROR SUMMARY." 

When Module Error Summary diagnostics appear in conjunction with 
coding errors, they can usually be ignored until all coding errors have 
been corrected. This is due to the fact that the Module Error Summary 
diagnostics were generated as a result of the incomplete processing of 
the statement containing the coding error. 

A complete list of the MARSYAS diagnostics appears on the follow- 
ing pages. 

\ 
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V B. MARSYAS DIAGNOSTICS 


Statement Error Messages 


NUMBER OF LEFT AND RIGHT PARENTHESES DOES NOT 
MATCH. 

Any parentheses that are opened must be 
closed. 

IMPROPER PUNCTUATION. A COMMA IS EXPECTED. 

IMPROPER PUNCTUATION. A $ SIGN IS EXPECTED. 

The statement has not been completed 
properly. 

IMPROPER PUNCTUATION. A COMMA OR $ SIGN IS 
EXPECTED. 

IMPROPER PUNCTUATION. A LEFT PARENTHESIS IS 
EXPECTED. 

IMPROPER PUNCTUATION USED IN TERMINATING 
SEQUENCE. 

EXPRESSION CONTAINS TOO MANY RIGHT PARENTHESES. 

IMPROPER PUNCTUATION. A SEMICOLON IS EXPECTED. 

WORD EXCEEDS 36 -CHARACTER LENGTH LIMIT. 

WORD BEGINS WITH IMPROPER FIRST CHARACTER. 

A MARSYAS word must begin with one of the 
characters ABCDEFGHIJKLMNOPQRSTUVWXYZ0123 
456789 . Words which are not numbers must 
begin with an alphabet letter. 


WORD CONTAINS AN INVALID CHARACTER. 

WORD CONTAINS IMPROPER NUMBER OF CHARACTERS. 

OPERATOR WORD IS NOT A MARSYAS OPERATOR. 

OPERATOR SHOULD BE INPUTS (OUTPUTS). 

OPERATOR WORD MAY NOT BE USED MORE THAN ONCE IN 
A GIVEN MODULE. 

EXPECTED OPERATOR NOT CORRECT. 

FUNCTIONAL DATA BASE PASS WORD INCORRECT OR 
MISSING. 

EQUATIONS MUST HAVE TWO SIDES. 

INVALID OPERATOR IN PRESENT CONTEXT. 

ONLY VARIABLES (NOT TIME) MAY BE DIFFERENTIATED. 

MNEMONIC WORD IS NOT A MARSYAS MNEMONIC; 

INTEGRATION MNEMONIC (CONSTRAINT) IS NOT A 
MARSYAS MNEMONIC (CONSTRAINT). 

MNEMONIC IS NOT A MARSYAS ARITHMETIC RELATION. 
The MARSYAS arithmetic relational operators 
are EQ, NE, GT, LT, GE, and LE. 




V B. MARSYAS DIAGNOSTICS 


Statement Error Mes sages 


WORD IS NOT A MARSYAS LOGICAL OPERATOR. 

The logical operators are AND, OR, and NOT. 

IMPROPER GRID TYPE MNEMONIC. 

IMPROPER SOLUTION MNEMONIC 

IMPROPER SAMPLE MNEMONIC. 

IMPROPER NUMBER OF PARAMETERS GIVEN FOR 
ELEMENT. 

IMPROPER PUNCTUATION. A COLON IS EXPECTED. 

IMPROPER NAME TYPE USED AS ARGUMENT. 

The statement requires a different type of 
name in the argument position. 

ARGUMENT SHOULD NOT BE DEFINED IN NAMING 
STATEMENT. 

The "old name" in a NAMING statement argu- 
ment list and the "new element" in a 
SUBSTITUTE statement should not be names 
defined in a NAMING statement. 

THE ARGUMENT HAS BEEN USED PREVIOUSLY IN THIS 
OPERATOR CONTEXT. 

Certain operations may be performed on a 
given argument only once. 


TIMESTEP CONSTRAINT SHOULD BE SPECIFIED. 

If the Runga-Rutta integration procedure is 
used, an integration step size must be 
specified. 

ORDER OF A TRANSFER FUNCTION MAY NOT BE VARIED. 

MORE THAN ONE WRITE PARAMETER HAS BEEN SPECI- 
FIED IN THE MODULE. 

A POST-PROCESSING module may write output 
at only one sample interval. 

CANNOT CHANGE GAINS WHEN THERE ARE LINEAR 
LOOPS. 

IMPROPER NUMBER OF INITIAL CONDITIONS. 

IMPROPER ELEMENT TYPE IN SUBSTITUTE STATEMENT. 

EXPECTED ARGUMENT IS INCORRECT TYPE. 

IMPROPER NUMBER OF ARGUMENTS. 

NAME HAS BEEN DEFINED MORE THAN ONCE. 

A MARSYAS name should be defined in only 
one statement. 

NAME IS USED INCONSISTENTLY. 

Characteristics which are implied by previ- 
ous uses of the name are inconsistent with 
the present use of the name. 




142 


t 


V B. MARSYAS DIAGNOSTICS 


Statement Error Messages 

PARAMETER HAS BEEN GIVEN TWO VALUES. 

CONTINUATION NUMBER REQUESTED NOT ON CONTINU- 
ATION FILE. 

NAME IS NOT IN MODEL DICTIONARY. 

END STATEMENT SENTINEL MUST MATCH BEGIN DECK- 

MODEL NAME ALREADY EXISTS IN FUNCTIONAL DATA 

NAME. 

BASE. 

A TERMINAL OF AN INCLUDED MODEL HAS NOT BEEN 

MODEL DOES NOT EXIST IN SYSTEM. 

CONNECTED. 

SOURCE CODE STATEMENTS ARE MISSING. 

FUNCTION (MATH REFERENCE) HAS INCORRECT NUMBER " • 

The program is incomplete. 

OF ARGUMENTS. 

SOURCE PROGRAM CONTAINS TOO MANY STATEMENTS. 


A MARSYAS deck may contain 99,999 state- 

USER DEFINED FUNCTION NOT ALLOWED AS A MATH 

ments. 

REFERENCE. 

THE SIMULATED SYSTEM IS TOO LARGE. 

NUMERIC FIELD CONTAINS AN IMPROPER CHARACTER. 

The MARSYAS system allows for 2-^-1 distinct 

Besides the digits, a numeric field may con- 

sequence numbers (objects). 

tain a sign, "E" or and only these 

characters. 

MODEL IS NOT INCLUDED IN THE SIMULATED SYSTEM. 

THE VALUE OF THE NUMBER IS TOO LARGE OR TOO 

FIRST STATEMENT OF PROGRAM SHOULD BE A 

SMALL. 

MARSYAS BEGIN. 

IMPROPER USE OF NUMERIC "E" FORMAT. 

PROGRAM NAME IS INCORRECT. 

THE NUMBER USED IS NOT PERMISSIBLE. 

MODEL INCORRECTLY RECORDED IN DESCRIPTION 

Certain numbers must be integers (non- 

TEMPORARY FILE. 

negative, positive, non-negative integers, 
etc.). Element subscripts and parameter 

FUNCTION WORD IS NOT A MARSYAS FUNCTION. 

numbers should be positive integers less 
than 2^-^. 


i 
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Statement Error Messages 

LENGTH OF PRINTED LINE EXCEEDS THE LINE LIMIT. 

SYSTEM HAS NOT BEEN PREVIOUSLY SIMULATED. 

An output line may be at most 132 characters 


long. 

SYSTEM HAS BEEN PREVIOUSLY SIMULATED OR 


CONTINUED. 

NUMBER IMPROPERLY FORMATTED. 



THE CONSTRAINT ILLEGAL IN PRESENT MODE. 

A TRANSFER FUNCTION PARAMETER B(N) IS ZERO. 


MODULE BEGINS WITH IMPROPER OPERATOR. 

- 

THE OPERATOR IS NOT PERMISSIBLE IN THE PRESENT 


MODULE. 


AN END STATEMENT MUST PRECEDE A NEW MODULE. 


A module in a MARSYAS program is completed 


with an END statement. A new module may 


begin after the END statement. 


THE SEQUENCE OF MODULES INVOLVING THE PRESENT 


MODULE IS INVALID. 


The order of modules for a given model 


must be chronological, e.g., a model must 


be described before it is modified or 


simulated. 


THE MODULE CONTAINS AN IMPROPER COMBINATION OF 


IF STATEMENTS. 
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Module Error Messages 

IMPLIED INTEGRATION AND DEFAULT INTEGRATION 

THE FOLLOWING CONNECT STATEMENTS CONTAIN 

METHODS DO NOT AGREE. 

IMPROPER OR AMBIGUOUS CONNECTIONS : 

MINIMUM TIME STEP LESS THAN STEP SIZE. 

THE FOLLOWING TERMINALS HAVE BEEN CONNECTED AN 
IMPROPER NUMBER OF TIMES; 

NO EXCITATION FUNCTIONS HAVE BEEN SPECIFIED 
FOR THE FOLLOWING INPUT TERMINALS; 

THE INITIAL CONDITIONS FOR THE FOLLOWING TRANS- 
FER FUNCTIONS ARE INVALID: 

EXCEEDING SYSTEM CONSTRAINTS CREATED AN ABORT 
SITUATION. 

- 

NO MODEL OUTPUT TERMINALS HAVE BEEN SPECIFIED ' 
FOR THE MODEL. 


NO ELEMENTS HAVE BEEN SPECIFIED FOR THE MODEL. 


THE MODULE CONTAINS NO PRINT OR PLOT STATEMENTS. 


A TRANSFER FUNCTION ORDER HAS BEEN CHANGED 
ILLEGALLY. 


THE FOLLOWING MARSYAS NAMES HAVE NOT BEEN PRO- 
PERLY DEFINED IN THE ABOVE MODULE: 
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EXAMPLES OF THE USE OF MARSYAS 
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VI. EXAMPLES OF THE USE OF MARSYAS 

A. Mechanical Extension Device - Example A 

An analysis of the motion for a mechanical extension device has 
resulted in the mathematical model presented in Figure 15. The motion 
of the device is represented by the quantities , # 2 > 6 3 an ^ 

which are the required outputs. The system input, $ , is subjected to a 
ramp input with slope 25.136 units. 

The system differential equation, as shown in Figure 15, is of 
the matrix-vector form 

a( e ) e + b( e x e ) 2 + c e + du - o (i) 


This system of equations is solved for 0 as 


e 


-1 


( e ) 


b( e x 9 r - c e 


DU 


] 


( 2 ) 


There are two methods of solving Equation (2) using MARSYAS: the 

block diagram method and the EQUATION method. The block diagram solu- 
tion of Equation (2), shown in Figure 16, makes use of specially con- 
structed DEVICES to perform the required matrix multiplications. The 
EQUATION approach solves a collection of scalar differential equations 
in conjunction with one DEVICE to perform the multiplication by A' 1 . 

The computer printout of the MARSYAS programs using each method is shown 
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VI Ao MECHANICAL EXTENSION DEVICE - EXAMPLE A 



WHERE 


(°-'l ) 
b cos (5 
c cos (03-0| ) 
d cos (04 -0, ) 




c cos (O 2 -O 2 ) 
d cos (0 4 -Q 2 ) 


- 0 


c cos (0 3 -0]) 
c cos (0 3 -02 ) 
(c- i 3 ) 

d cos (04 -0 3 ) 


d cos (0 4 - 0| ) 
d cos (04-02 ) 
d cos ( 04-03 ) 

(d-i 4 ) 


b sin (0 2 -0] ) 
c sin (0 3 -0 1 ) 
d sin ( 04 - 0 ! ) 


-b sin (0 2 - 0, ) 

0 

c sin (0 3 -0 3 ) 
d sin (04 _0 2 ) 


-c sin (0 3 -0] ) 
-c sin ( 0 3 - 02 ) 
0 

d sin 


-d sin (0 4 -0]) 1 
-d sin ( 04 „ 0 2 ) 
-d sin (0 4 _0 3 ) 

0 


C= - 


I 

l 2 


2k, -k, 
-k 0 2k 0 


0 *>k 3 2k 3 _k 3 

0 0 -k4 k 4 


D= — - 


l 2 


-k, 0 

0 0 

0 0 

0 -k 4 


FIGURE 15 MATHEMATICAL MODEL OF MECHANICAL EXTENSION DEVICE 
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VI A. MECHANICAL EXTENSION DEVICE - EXAMPLE A 


on the following pages. Note that the translation time for each module 
is automatically printed out with the coding. 



VI A. MECHANICAL EXTENSION DEVICE - EXAMPLE A 
(Block Diagram Option) 

a t g i n : £ x a fi p l e a « 

I) t V I C t : A I N V f 3 » 4 > 7 s 

SUBROUTINE A I N v { IJ , Y , P ) 
oimn'isiUN U(3) s r(4),p(7) 

[) I MtNS j UN A I 4 » 4 ) , v ( 3 ) s J C ( 4 ) .* A ( 4 ) 

A ( 1 s 1 1 = P ( 1 ) 

A(!»2)=tlb)oCUS(U(2)”Ull)l 
A I 1 , 3 ) = P ( 6 I o C 0 S ( U l 3 ) - 0 ( 1 ) > 

A ( l » 4 ) = H ( 7 ) » C 0 S ( U ( 4 ) » U < 1 ) I 
A I 2 » 1 > = A ( 1 , 2 ) 

M2, 2 > = P I 2 ) _ . 

A ( 2 , 3 ) = P ( 6 ) o 0 OS ( U ( 3 > -U ( 2 ) » 

A ( 2 » 4 ) = P ( 7 ) * C 0 S < IJ I 4 ) “ U I 2 ) ) 

A ( 3 » 1 ) » A ( l , 3 ) 

A(3»2)=A12,3) 

A ( 3 » 3 > = P ( 3 > 

A ( 3 ,4) « P ( 7 ) * C 0 SJ U (J ) ~ IJ l ( 3 ) J 

Alt, l l=A( Ml 

A ( 4 , 2 ) = A I 2 , ‘I ) _ 

A ( 3 , 3 ) s A ( 3 , 4 ) 

A ( 4 , 4 I = P ( 3 ) 
is) = 3 

V (12 = 3.0 

CALL G J* I A , N , N ,N , N , S 1 5C , JC , V ) 

DO 2U0 1 = 1 »4 
2 0 0 X C I > ~ U ( 4 ♦ 1 ) 

M a 4 
Na 1 

C A L L _M X M l. I I A , X , r , M_, M , M , M , M ) 

G 0 1 0 4 0 U 


1 BO 

IP I V i 1 

J ) 132,151 ,152 

1 B l 

UP 1 IP- (6,133) 

1 3 3 

P o w m a r 

I * MATRIX A I S 


SIDP 


1 3 2 

WHITE ( 

6,154) 

I 3 4 

format 

I * OVEHFLO'V * ) 


STOP 


4 00 

K £ T U K N 



E No 

DEVICE: HD, 8 » 4 , 35 

S VdKOUTjNt R B ( U , Y , p ) 

0 1 ME N S I U N U ( 8 I , Y I 4 ) , P ( 3 ) 

D 1 MENS I DM [3 ( 4 , 4 J , X ( 4 ) 

SI I , I 1=0, 

B I 1 t 2 I = P ( 1 ) * S 1 N ( D ( 2 ) - U ( t 1 ) 
B ( l , i ) = P(2)«S!N(U(3]-U(l I ) 
__B ( I » H > a P I 3 ) « S l N ( U (Jt ) -U I 1 M 
B'l 2 , l » =-B ( V , 2 > 

B ( 2 »2 I =0° 

B ( 2 t 3 ) = P I 2 ) « 5 J M l U ( 3 ) ~ (J I 2 1 ) 
I) ( 2 , 4 I = P | 3 ) c S I D ( U ( 4 ) ~ U ( 2 J ) 
B I 3 , 1 ) a - □ I 1 , 3 ) 

_B l 3 , 2 ) a -H ( 2 a 3 ) 

B ( 3 * 3 ) = rj » 

B ( 3 , H ) = P I 3 ) « S I N ( U I 4 ) - U ( 3 ! > 
O I 4 , l ) ~ - B i l , 4 ) 


150 



VI A. MECHANICAL EXTENSION DEVICE - EXAMPLE A 
(Block Diagram Option) 

Hi '4 , 2 ) = - 0 ( ?. , H ) 

►M *4 V *U = o e 

DO IUU 1 = 1,*) 

1 U 0 X ( I ) " U ( 1 + ‘3 ) a o 2 
Ms') 

Ns 1 

CALL M X •11.. 1 ( H , X , Y , N , M , N , M , rl ) 

RtrUKH 

t N L) 

_ DEV ICE; _ M C , ■), H , HSi 
suh ko u r j tit mc ( u , y (Pi 

D 1 ME US I OlO U ( M ) , Y ( 4 ) , P ( 4 > 

K E A L K 

I) I M E N 5 1 U N K ( 'i , '( ) 

K ( 1 * n =2«P ( 1 ) 

K ( l , 2 ) ~-P ( 1 ) _ _ 

k i r. 3 ) =o« 

K ( 1 » H ) « 0 s 
K ( 2 * 1 ) =~P ( 2 ) 

K(2,2)=2oP(2) 

IU3t3)=-P(2) 

__K ( 2 . H ) s 0 • 

K ( 3 t 1 } ™ (J b 

Kt3,2)a-P<3> _ 

K (3,3)=3«P|3) 

K ( 3 , <4 ) s - P ( 3 > 

K 1 H i 1 l-Oi 
K ( H i 2 ) = 0 « 

K"( ^ , 3 ) =~P ( j~ 

KlHitlsPco 
M s H 
Ns l 

CALL IXHLT (K,U,Y 
KEfUKN 

‘ END . 

device: OH, 2 * M , 2 s_ 

‘ SUbKOUT t IVE DM ( U , Y , P ) 

O I II £ N S l O N U ( 2 > , Y { h ) , P ( ^ * 

DIMENSION D ! ■) i ? ) 

D ( 1 , 1 ) = - P ( I ) 
p -■ 2 j * 0 • • 

D ( 2 « t ) =0. 

" ■ D I 2 t 2 J = EJ» 

D ( 3 » 1 ) ~U« 

„ ( y;^ ) =u » ~ “ “* 

D ( s « i ) = 0 • 

= - P I 2 ) — — 

M = 'l 
N a l 
t.a2 

CALL M X ML T { 0 , U , Y , M , L. , N , M , L ) 

R E T U K N 
E N D “ 

f-1 0 D E L : E X T E N S I 0 M DEVICES 

INPUTS: THE I AO, DELTAS 
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Output s : 
t LLrie.N r s : 


C 0 M |\| t C ! 


£ N D <6 


S l M U L A r E : 

EXC 1 Tfc : _ 

sToPi f- : 

E N [.) *5 

ALL I N I r I At C 


p h in r : 
SAMPLE : 
END'S 


end : 


VI A. MECHANICAL EXTENSION DEVICE - EXAMPLE A 
(Block Diagram Option) 


theta 
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T A 

3 
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the 
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3 
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A l N V , 


A X ( 
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fe 
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... 
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o 

) s 















AD , 


AO l 

> 

AO 2 , 


A 

0 3 

9 

/A 

03s 












IN, 


I N I 

9 

i :i2 , 


I 

*13 

9 

1 


I 

N 

5 9 

1 

H 

6 

9 

1 

N 

7 4 

I N 8 si 

theta 

u 

, I •> 0 

X 

tt I , A 

0 

1 . 

9 . 

b 

tt A 

X* 1 

9 . 

I 

N 1 

9 

i 

N 

2 

9 

T 

H f - 

T A ). S . .. .. 

pelt a 

* 

2 VOX 

» 

2 » AO 

2 

9 

6 


AX 

tt 2 8 

I 

N 

3 9 

I 

N 

8 

9 

T 

H 

ET. 

A? S 

. OX tt 3 , 


A 0 3 s 

7. 

tt A x a 3 

V . 


.1 >' 

3 

9 - 

IN 6 

?„ 

T 

HE. 

T A 

3 

<5 





. 

DXH4 , 


AUH, 

8 

tt A X tt 3 

9 


I ft 

7 

9 

I N g 

9 

T 

HE 

TA 

3 

5 






l N 2 , 

1 

u c a a I 

9 

AD 1 £ 
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_ 


.. 

. _ 


_ _ _ . _ 

I N 3 t 

2 

ti C X » 2 

• 

A 0 2 s 


















IN 6 , 

3 

rt C X * 3 

9 

A D 3 S 


















I NO , 

3 

S C X H <4 

6 

A D 3 3> 


















JJT 2 s ~ 

1 

SAX V _ 




_ _ 





... 

_ _ . 







. . — . 






_ 




_ „ 

I N 8 , 

2 

» A X 5 




















I N 6 , 

3 

ft A X S 




. 




.. 








_ 


. 


.. . _ . . 

I MO i 

3 

it A X * 




















.IN?., 

I 

it 13 X « I 

9 

Apis 


















I N4 , 

2 

tt B X * 2 

9 

A 0 2 S 


















I N 6 4 

3 

tt B X Jt 3 

9 

A 0 3 S 


















f N 0 4 

*♦ 

h IJ X « 3 

9 

ACT® 


















INI, 

S 

» 3 X <S 




















i N 3 , 

6 

PSX"' 




















I NS 4 

7 

ft 3 X * 




















I N 7 , 

0 

8BX<5 





















MODEL l ♦ 5 l 3 3 G 0 SEC. 

EXTENSION DEVICE* 

J)ELTa » FK aMP^S* 1 36 ) , THEJ aOj ZEROS _ 

T l M E • <3 T ° 6 2 5 S ■ 

0 N 0 I T i U N S HAVE BEEN LIE K A U L I £ D T 0 7.E RO„ 

SIMULATE To VM'6.00 SEC", 

D £ L I' A t THE TaJ. i r tL K T *1? »_T HE T A 3 JHE TA 3 S 

step, is ” "" ~ 


POST PKOC _ 3. 4 3 3 2 30 SEC_„ 

TxaDPLEAS 

total lTTvsToo " 'sec . 
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VI A. MECHANICAL EXTENSION DEVICE - EXAMPLE A 
(Equation Option) 


B EG 

i m : 

■ - 

-- • - 

EXAMPLE A 

«r> -- 


.... 

DEV 

ICE 

D 


AINv, B c 



7 $ 

S 

IJ R R 

nu 

TIME 

A I N V ( U , Y , 

R ) 



D 

I ME 

us 

I'M u 

( S) » v ( H i , 

P ( 7 

) 


n 

I M E 

NS 

ION A 

( 4 , ’1 ) t V ( 3 

) , U 

c 

( M ) , X 1 M) 

A 

( 1 » 

1 ) 

RP ( M 





' "A 

( 1 . 

2 1 

= P< 5 ! 

o C 0 S { IJ ( 2 1 

-U( 

; 

) ) ' 

A 

< 1 » 

3 ) 

= P ( 6 ) 

* C 0 S ( U ( 3 ) 

U ( 

1 

) ) 

A 

< 1 . 

M ) 

= P ( 7 ) 

o c o s ( u ( q i 

-U ( 

I 

) ) “ ‘ ' ~ 

A 

( 2 » 

1 ) 

= A ( 1 , 

2 ) 




A 

( ? , 

2 ) 

* P ( 2 ) 




“ ’ 

A 

( 2 , 

3 ) 

= P ( (> ) 

c C 0 S ( U ( 3 ) 

~U ( 

2 

) ) 

A 

(? , 

M ) 

8 P ( 7 ) 

©cos (u ( m ) 

~ u ( 

2 

yi 

A 

< 3 , 

1 ) 

*A ( 1 , 

3 ) 





A ( 3 , ? ) = A ( ? , 3 ) 

A ( 3 , 3 ) = P ( 3 ! 

A ( 3 9 M > = P ( 7 ) © C 0 S ( IJ { M ) - !,! ( 3 J ) " ' ' 

A ( M , 1 ) = A ( 1 , /( ) 

A ( 'Vi ? ) s A ( 2 , <» ) 

A(M,3)=A(3,q) 

A ( M , m ) - P ( M ) 

n = m 

v u ) = 3 o c; 

CALL G JK ( A , N , N , N , N , <$ ! GO , JC , V ) 

. . — D o- 2 no 1 = 1 , q : 

Z qo x ( i ) * u ! m * i ) 

11*1 
H* 1 

CALL MX'iLTf 
GO TO ^"O 

■ ISO' IE (VMM 1 5 2 * 1 ' *3 1 ,15 2 

151 WRITE (6,153) 

|15 3 FORMAT ( • MATRIX A IS SINGULAR* ) 

STOP 

| 5 2 WRITE (6,1 SR) 

} 5 M FORMAT ( * OVERFLOW*) 

MOO RETURN 
E N n 

model: extens i omdev t cf, , equations 

inputs: to, dels 

outputs: t o i , top, tot, toms 


• . . — . 

• 

' - PT 

■8--B 

«■ SI 
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J © • 
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- c — 1 -<> — 2“+' 
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T 1 * *> e 2 
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( 1 Q 
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VI A. MECHANICAL EXTENSION DEVICE - EXAMPLE A 

(Equation Option) • — - — - 

♦' D a. A o S , 

5 K a 1 ,3 O <5 

; TO t = T1 -5 - - - - -- - - 

5 T O 2 = T ? <5 

: T 0 3 = T3 S - - - - - - - - 

: T o '{ s T -4 <3 

-T t “• * *-Af MVt A,«I233«S- , A o 3 A 1 ?7S , A 0 ?99 1 AS V~ 3o?92CS5 , 

A e 3 A A 2 7 S , A a 3C 9 1 AS , A e 2 9 7 OSS , T ! , T2 t T 3 , 

T*» , PI ■ P2 ' , P3 , P 9 , 1 ) s - — 

: T? 6S = A 1 1 J V ( A 0 2 3 3 E - 5 » Ao.3A127S , A . ? 9 9 I A S , 3 «, 292Q55 , 

• A.3AA?75 , A»3091AS » Ao29?e5S , T 1 T 2 , T3 v 

Tq,pi,p?,P3,p*»,?)s 

T 3 - ♦-* " «r*# T N V C A ^ « 2 3 3^ 5" 6' « 3 A !' 2 7 5 . ' A i 2 9 9 1 A 5 ' * * 3 <> 2 9 2 0 S 5 ' , 

Ao3AA?7S » Ac 309 IAS , 6 <,29 2055 , T) , T2 • T3 , 

" ' ■■■“■ T'4 , Pi , -P2 , P 3 , P9 ,• "3 ) S — 

: T9 = AINV(A.*J?33R5 , Ao 3A!275 » A 0 ?99l6S , 3o2920S5 , 

A o 3 A A 7 7 S , A.3CM1AS , Ao292H55 , ”TI "» T2 . '~'T3 V' 

^ _ _ T1 , P| , P2 , PJ , pi , !| ) , 

SIMULATE! EXTENSION DEVICES 


t x C I T E : DEL,, PR AMP ( 2B o 1 3A ) S TO , ZEROl 

s T 0 p J F ; TIrlF„GTe»?B$ 

ALL INITIAL CONDITIONS HAVE BEEN DEFAULTED TO ZERO, 
SIMULATE M»3CA9E'0 SEC, 

print: del, rot , to?, -tot, to'h-- - 

Sample: step, is 

post proc . . . 9 P 3 c > 8 9 0 o sec„ 

end: e x a h p l e a * 

TOTAL I 1 , 9B 1 800 SEC c 
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VI B. DO NOTHING SYSTEM WITH SUBMODELS - EXAMPLE B 


✓ 


B. Do Nothing System with Submodels - Example B- 

The system block diagrammed in Figure 17 contains most of the ele- 
ments contained within the MARSYAS library. One purpose of presenting 
this fictitious model is to illustrate the use of submodels and how 
they are interconnected within the main model. Most of the features of 
MARSYAS described on the preceding pages (with the exception of EQUATION) 
have been incorporated in this example. It is intended to be used as a 
learning aid for the beginning MARSYAS user. 
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BETA 




























VI B, 


MODEL OF DO NOTHING SYSTEM 


BEGIN: EXAMPLE B$ 

* A USER DEFINED ELEMENT$ 

DEVICE: AT, 1, 1, 1$ 

SUBROUTINE AT(AI,AO) 

DIMENSION AI(1) , AO(l) 

AO(I) = ATAN (AI(1) ) 

RETURN 

END 

* A USER DEFINED FUNCTION$ 

FUNCTION, DFSIN, 5$ 

FUNCTION DFSIN (Dl, D2 , D3 , D4, TIMEX) 

C THIS FUNCTION IS ALSO PART OF THE MARSYAS LIBRARY 

IF (TIMEX .LT. D4) GO TO 1 
DFSIN = Dl * SIN (D2 * (TIMEX - D4) + D3) 

RETURN 

1 DFSIN =0.0 
RETURN 
END 


MODEL: DO 

INPUTS : 
OUTPUTS : 
ELEMENTS : 


NOTHING M0DEL$ 

11,12,13,14,15,16,17$ 

01,02,03,04,05$ 

AD, AD1,AD2,AD3,AD9,AD12$ 

SI, NL7 (1. ,2. ) , SI2(2. 0,3.0)$ 

PF, BL7 (1.0)$ 

ML, NLl$ 

RE, NL3$ 

TF, BL3(2,1. ,2. ,3. ,4. ,5. ,6.)$ 

HS, NL4(5,-2. ,-2. ,+2. ,-2. ,46. ,+2. ,0. ,+2. , +2., -2.)$ 

HS, HL4(5, -1.1, -2.0, +3., -2.0, +7'. ,+3,0. 0,0.0, 0. 0,0.0)$ 
CO, NL2(1. 0,2.0)$ : C02 (2 . ,3 . ) $ 

LM, PL4(2. ,-l. ,3.)$ 

DV, DV1$ 

IT, IT2$ 

CM, CM8 (3 . ) ,CM9 (4. ) ,CM10 (8. ) ,CM11 (5 . ) ,CM12 (7 . ) ,CM13 (6 , ) , 
CM15 (9. ).$ 

AT, ATT (1.3)$ 


SUBMODEL 

SUBMODEL 

SUBMODEL 

CONNECT: 


BETA; INPUTS: IB1,IB2,IB3; OUTPUTS: OBl,OB2,OB3$ 

GAMMA; INPUTS: IG1,IG2; OUTPUTS: 0G1,0G2$ 

ALPHA; INPUTS: IA1; OUTPUTS: OA2$ 

11, I#NL1, ADI, NL4, AD2, PL4 , IGl, OG2 , AD3 , BL7 , NL7 , 01$ 

12, 2#NLl$ : NL4, CM10, AD1$ : NL7 , CM8 , 1#DV1, AD3$ 

BL7, CM9, 2#DVl$ 
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VI B. MODEL OF DO NOTHING SYSTEM 


CONNECT: 14, 1#NL3#1, ADl$ : 16 , 3#NL3$ 

: 15, 2#NL3#2, IA1, OA1, IG2$ : OA1, IB3$ 

: OA1 , 1# IT 1 , ATT, 05$ 

: ATT, AD9$ : IT2, CM12 , AD9$ 

: AD9 , SI2, NL4, CM13, AD9$ 

: NL4, C02 , CM11, 2#IT2$ : OGl, 04$ 

CONNECT: 17, IB2$ : 13, AD12, NL2 , CM15 , AD12$ 

: NL2 , AD2$ : OBI, 02$ 

: NL2 , BL3 , IBl, 0B2 , AD3$ 

: 0B3 , 03$ 

END$ 

MODEL: GAMMA$ 

INPUTS: IG2 , IGI$ 

OUTPUTS: 0G2 ,0G1$ 

ELEMENTS: IT, NL6$ : AD, AD5$ : TF, BL6(4, 0.0, 1.0, 0.0, 0.0, 

1.0, 1.0, 3.0, 2.0, 1.0)$ 

: CA, CM1 (2 , +1)$ 

CONNECT: IG1, 1#NL6, AD5 , BL6 , 0G2$ : BL6 , CMl, AD5$ 

: NL6, 0G1$ : IG2 , 2#NL6$ 

END$ 

MODEL: BETA$ 

INPUTS: IBl, IB2 , IB3$ 

OUTPUTS: 0B3 , 0B2 , OBl$ 

ELEMENTS: AD, AD4$ : DV, NL5$, TF, BL4(4, 0.0, 1.0, 0.0, 0.0, 1.0, 1.0, 

0.0, 3.0, 0.0, 2.0)$ 

: AT, BL5$ 

SUBMODEL: GAMMA; INPUTS: IGl, IG2; OUTPUTS: OGl, OG2$ 

CONNECT: IBl, AD4, 2#NL5 , BL4, IGl, OG2 , OB2$ : BL4 , BL5 , AD4$ 

: IB2 , 1#NL5$ : NL5 , OBl$ : IB3, IG2$ : OGl, OB3$ 

END$ 

MODEL: ALPHA$ 

INPUTS: IA1$ 

OUTPUTS: 0A1$ 

ELEMENTS: CM, CM2$ : DS, BL2(2.,-3.)$ 

CONNECT: IA1, CM2, BL2 , 0A1$ 

END$ 

SIMULATE: DO NOTHING MODEL$ 

EXCITE: 17, DSTEP (1.0, 2.0)$ 

: II, ZER0$ : 12, FSTEP(1.0)$ : 13, FSIN(0. 8,60.0,0.5)$ 

: 14, DSTEP(2. 0,5.0)$ : 15, FRAMP(2.0)$ 

: 16, DAMP(1. 0,2. 0,3.0)$ : 17, FPULSE(1. 0,2. 1)$ 

INTEGRATE: RK, TIMESTEP, 0.01$ 
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VI B. MODEL OF DO NOTHING SYSTEM 


CHANGE: RUN2 ; EXCITE: 17, DFSIN(1.0, 1. 0,0. 0, 1. 0) $ 

CHANGE: RUN3 ; PARAMETERS, GAMMA, BL6(3, 2.6)$ 

TERMINATE IF: TIME .GE. 20.0$ 

END$ 

PLOT: LOG LOG (0.0,20.0,1) , 01, 02, 03, 04, 05$ 

PRINT: II, 12, 13, 14, 15, 16, 17$ 

SAMPLE: STEP, 10$ 

FOURIER: 04(8.0)$ : 06(8.7)$ 

END$ 

END: EXAMPLE B$ 



VI C. VEHICLE STABILIZATION SYSTEM - EXAMPLE C 

C. Vehicle Stabilization System - Example C 

The Vehicle Stabilization System shown in Figure 18 is an example 
of a system containing nested submodels and differential equations. 
This is an example of mixing the block diagram and Equation methods to 
solve a problem. 



VIC. VEHICLE STABILIZATION SYSTEM-EXAMPLE C 



CONTROL SYSTEM XI 




FIGURE. 18 VEHICLE STABILIZATION SYSTEM BLOCK DIAGRAM 


< A 












VI C. VEHICLE STABILIZATION SYSTEM - EXAMPLE C 


BEGIN: VEHICLE STABILIZATION SYSTEM$ 

MODEL: CONTROL SYSTEM XL$ 

INPUTS: INI, IN2$ OUTPUTS: HI, VI $ 

ELEMENTS: RE, RESOLVER $ 

SUBMODEL: ACTUATOR STAGE 1 AND 3; INPUTS: Al, A2, A3; 

OUTPUTS: ACTI, ACT2$ 

SUBMODEL: GIMBAL3; INPUTS: Gl; OUTPUTS: GIML, GIM2$ 

CONNECT: INI, Al, ACTI, l#RESOLVER#l , Hl$ 

: ACTI, 2#RESOLVER#2 , Gl, GIM1, A3$ 

: IN2, A2, ACT2 , 3#RESOLVER $ 

: GIM2, VI $ 

END$ 

* THE FOLLOWING IS THE DESCRIPTION MODULE FOR SUBMODEL ACTUATOR STAGE 1 AND 3$ 

MODEL: ACTUATOR STAGE 1 AND 3$ 

INPUTS: Al, A2 , A3$ OUTPUTS: ACTI, ACT2$ 

ELEMENTS: AD, ADI, AD2$ : TF, MOTORA(2, 0.0, 3.0, 5.0, 2.0, 4.0, 7.0)$ 

: LM, LIMITERC(1 . , -3., 3.)$ 

SUBMODEL: GIMBAL3; INPUTS: Gl; OUTPUTS: GIML, GIM2$ 

CONNECT: Al, AD2$ : A2, ADI, MOTORA, AD2 , Gl, GIML, ACT1$ 

: GIM2, ACT2 $ : GIM2, LIMITERC, AD1$ : A3, AD1$ 

END$ 

* THE FOLLOWING IS THE DESCRIPTION MODULE FOR SUBMODEL GIMBAL3$ 

MODEL: GIMBAL3$ 

INPUTS: Gl $ OUTPUTS: GIML , GIM2$ 

ELEMENTS: AD, AD1$ : CM, GAIN5 (-10.3) $ 

SUBMODEL: CONTROLLER; INPUTS: X ; OUTPUTS: Y, Z$ 

CONNECT: Gl, ADI, X, Z, GAIN5 , AD1$ : Z, GIM2$ : Y, GIML $ 

_ END$ 

* THE FOLLOWING IS THE DESCRIPTION MODULE FOR SUBMODEL CONTROLLER$ 

MODEL: CONTROLLER, EQUATION$ 

INPUTS: X$ OUTPUTS: Y, Z$ 

EQUATION: Z' + Z = X$ 

: 5.0 * Y" + 2.0 *Y' + 3.0 * Y = 2.0 * X - 5.0 * X ** 2 

END$ 

SIMULATE: CONTROL SYSTEM XI $ 

INITIALIZE: ACTUATOR STAGE 1 AND 3, MOTORA(1.5, 12.0)$ 

EXCITE: INI, FSTEP(5.0)$ : IN2, FSIN (1.0, 3000.0, 0.0)$ 

TERMINATE IF: TIME .GT. 2.0$ 

END$ 

PRINT: (VEHICLE STABILIZATION SYSTEM, RUN ONE) $ 

: INI, IN2, HI, VI $ 

FOURIER: HI (1.3)$ : VI (1.3)$ 

END$ 

END: VEHICLE STABILIZATION SYSTEM$ 
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o(f) - 


n iS THE NUMBER OF INPUTS 
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CLASS 

BLOCK DIAGRAM 
SYMBOL 

# OF 
INPUTS 

# OF 

OUTPUTS 

MNE- 

MONIC 

1 ! 

INPUT-OUTPUT LIST OF PARAMETERS IN THE ORDER IN WHICH j 

RELATION | THEY APPEAR IN THE ELEMENTS STATEMENT j 

1 1 

EXPG- 

NENTIATOR 

- 

EXP 

o(t) 

1 

1 

EX 

II 

o 

' ' Jl 1 

i 

i 

mnE J 

1 

LOGARITHMIC 

i(t) _ 

4^-' 

LOG 

^J> 

1 

1 

LN 

o(f) = log e i(t) , i(t) > 0 

1 

I 

NONE 

SINE 

i(t) _ 

SINE 

o(tL 

1 

1 

SI 

o(t) = sin(ai(t) + b) 
(^assumed to be in radians) 

(a,b) 

i 

i 

i 

. 

COSINE 

i(t) 

-—{> 

COS 

«>(*) 

1 

1 

CO 

c(t) = cos(ai(t) + b) 

(]_ assumed to be in radians) 

(0,b) I 

1 

1 

i 

■ 


TAN 

o(t) 

E>- 

1 

1 

TN 



o(t) = tan(ai(t) + b) 
(^assumed to be in radians) 

. 

1 

i 

(a. b) j 

; i 

j 
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# OF 

OUTPUTS 

MNE- 
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INPUT-OUTPUT 

RELATION 
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THEY APPEAR IN THE ELEMENTS STATEMENT 

SWITCH 


o(t) 

2 

1 

sw 

i 2 (t), i T (t)>Tl 

o(t) - 

0, i|(i) < T 

T 

SIGN CHANGER 

T(t) oU)*" 

1 

1 

sc 

o(t) = - i(t) 

NONE 

TIME VARYING 
COEFFICIENT 


o(t) 

1 

1 

TV 

o(i) = i(t) °f(?) 

y 

(*2' y 2^ 

*> 

<*3 ■»:)>' 
MM'" 

(n, f 1 , y-j ,t 2 , y 2 ""tn , Y n ) 

1 

I 

1 
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INPUT-OUTPUT 

RELATION 


o(t) = I i (t)| 


o(t) = i(nT), 
hT5t<(n v 1) T, 
n = 0, 1, 2, ... 


o(t) = i i (t) • i 2 (*) 


o(t} = 


i,(t) 

i2( } ) 


IF SUCH A DEVICE IS INCLUDED IN THE SYSTEM, THE RUNGE KUTTA SCHEME 
FOR INTEGRATION MUST BE SPECIFIED IN AN INTEGRATE STATEMENT IN THE 
SIMULATION MODULE WITH TIMESTEP LESS THAN T„ 


LIST OF PARAMETERS IN THE ORDER !N V/KiCH 
THEY APPEAR IN THE ELEMENTS STATEMENT 
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DSTEr (d 1 , d 2 ) 


DFSiN (d j , d 2 ^ dj , 
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